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1.  GERESS  STATUS  REPORT 


January  1992  -  December  1992 


Michael  L.  Jost 


1.1  General 

The  German  Expenmentai  Seismic  System  (GERESS)  is  a  c(XH)erative 
research  program  of  Southern  Methodist  University  Dallas.  Texas,  and  Ruhr- 
University  Bochum,  Germany.  It  is  an  extensitm  of  the  Scandinavian  regional 
array  network,  i.  e.,  NORBSS,  ARCESS,  and  FINESA.  into  Central  Europe. 
The  GERESS  array  is  located  in  the  Bavarian  Forest  area  at  the  souiheasmm 
part  of  Germany  near  the  border  to  Austria  and  Czechoslovakia.  The  array 
consists  of  25  stadons  with  verticalocotnponent  short>period  Teledyne  Geotech 
GS-13  instruments  sampled  at  40  Hz.  In  addition,  four  of  the  sites  include  hor¬ 
izontal  component  instruments.  At  the  key  station  of  the  array,  GEC2 
(48.84511  N,  13.70156  £.  1132  m).  a  STS-2  is  sampled  with  80  Hz  (broad¬ 
band  element).  The  geometry  of  die  array  is  based  on  concentric  rings  provid¬ 
ing  an  overall  qiertore  of  about  4  km.  The  GERESS  ^lermre  is  larg^  than 
that  of  NORESS  by  a  facttv  of  4/3  which  resulted  fiom  signal  and  noise  corre- 
ladcm  measurements  by  Harjes  (1990).  The  array  became  fully  operational  in 
January  1991.  Data  fiom  the  array  were  cmuinuoosiy  transtxdtttd  to  NORSAR 
and  to  the  Institute  of  Geophysics.  Ruhr-Universi^  Bochum  via  64-kint  tele¬ 
phone  lines  in  1992.  In  early  1993,  the  NORSAR  line  has  been  terminated  ami 
a  wide  area  network  (WAN)  has  been  established  between  the  GERESS  data 
centers  in  Bochum  and  Kjeller  (Hgure  1-1). 


IJt  GERESS  Observatory  at  Rutar-University  Bochum 

The  Department  of  Geophysics  of  Ruhr-University  Bochum  operates  an 
experimental  on-line  processing  system  for  GERESS  data  (Harjes  et  aL,  1993). 
Tto  system  uses  software  developed  at  NORSAR  (Mykkeitveit  and  Bungum, 
1984;  Fyen.  1987).  During  the  Group  of  Scientific  Expens  Technical  Test 
(GSETT-2)  and  in  1992  numerous  updates  to  the  program  package  were  in^le- 
mented  in  Bochum. 

The  on-line  processing  consists  of  3  stqps:  detection,  flc-analysis.  and  loca¬ 
tion.  The  first  stage  of  the  on-line  processing  accesses  data  in  30  second 
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segments  and  runs  a  STA/LTA  detector.  The  detector  presently  recognizes  an 
onset  if  the  STA/LTA  ratio  for  a  filmxed  trial-beam  exceeds  a  threshold  of  4. 
The  next  step  of  the  on-line  processing  is  the  transfmmanon  of  a  3  second 
filtered  data  segment  at  each  onset  time  (derived  fixmi  the  ttetection  time)  into 
the  frequency-wavenumber  domain.  As  a  result,  the  slowness  and  l»u:k- 
azimuth  of  the  phase  is  determined,  from  the  slowness  infonnation,  seismic 
phases  are  identified.  The  final  step  of  the  mi-line  data  processing  is  the  loca¬ 
tion  of  events.  The  seismic  phases  as  identified  in  the  £k-analysis  are  associated 
to  events  in  this  step.  From  the  arrival  time  difference  of  regional  phases,  the 
distance  to  the  epicenter  can  be  determined  from  the  Jeffreys  >  BuUen  travel 
time  tables  for  regional  seismic  phases.  Together  with  a  mean  back-azimuth, 
the  qncenter  locations  of  local  and  regional  events  ate  determined.  Far  1992, 
18  events  have  been  automatically  located  each  day  on  average  (Table  1-1). 
Figure  1-2  shows  the  automatic  event  locations  for  1992. 

Since  October  1992,  analysts  started  to  review  regional  events.  They  lep- 
icked  onset  times,  redetermined  back-azimuths  and  phase  velocities,  and  relo¬ 
cated  events  using  the  LOCSAT  routine  which  is  a  modified  version  of  TTA- 
ZLOC  (Bratt  and  Bache,  1988).  Hgure  1-3  (Schulte-Theis,  pers.  comm.)  shows 
the  relocated  events  from  Oct-Dec  1992. 


Table  l-l 

Statistics  on  the  GERESS  Bochum  On-line  Processing  for  1992 
_ number  number/day 


detectkms 

62749 

171 

f-k  analyses 

57090 

156 

locations 

6540 

18 

At  Bochum,  local  magnitudes  are  calculated  by  die  following  procedure 
(the  algorithm  was  inq)lemented  by  J.  Schweitzer): 

1)  Deteraunation  the  optimum  S-phase  beam. 

2)  I£^-pass  filter  at  0.7  Hz. 

3)  Transfonnatioa  of  beam-trace  to  a  Wood-Anderson  record. 

4)  Pickiag  of  maximum  amplitude  of  the  S-waves  and  coirection  epicentral 
distance. 

The  ML  magnitude  distance  correction  from  Gutenberg  and  Richter  (1956)  was 
changed  following  the  procedure  of  Ebel  (1982)  using  die  attenuation  relaiton 
for  France  from  Nicolas  a  al.  (1982).  In  addition  to  the  beam  magninide,  local 
magnitudes  have  been  determined  for  the  single  station  C2  and  for  die  torizon- 
tai  componenis  in  a  similar  way.  Usually,  ^  3  magnitude  values  diffor  less 
than  aim  0.3  magnitude  units.  On  the  other  hand,  GERESS  local  magnitudes 
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appear  low  by  at  least  OJ  noagnitude  units  compared  to  values  published  by 
otho’  European  networks. 


Since  Sep.  1991.  the  results  of  the  on-line  processing  have  been  sent  to 
interested  institutions  via  e-mail  (e.g.  Bundesanstait  fur  Geowissenschaften  und 
Rohstoffe  (BGR),  Frankfun  University,  Graienberg  Array,  ORFEUS,  seismo- 
logical  institutes  in  CSFR,  Hungaria,  Slowenia,  SwitzeriW;  Stuttgart  Univer- 
si^).  In  addition,  the  On-line  Processing  Display  Manager  (lost,  1991)  has 
been  used  by  many  scientists  for  near  real  time  infonnatitHi  on  parameter  d^ta. 
Recently,  a  Data  Request  Manager  (DRM)  by  Krake  Inc.  has  been  installed  for 
trial  operadon.  The  DRM  should  enable  the  txansmissitm  of  data  (GSE,  SEED) 
and  interfaces  to  the  German  Regional  Network  of  Broad-Band  Stations. 

Since  late  October  1991,  alert  messages  of  strong  teleseismic  and  regional 
events  are  automadcally  sent  to  NEIC  (Golden,  CO)  and  BGR  in  near  real¬ 
time.  46  times,  GERESS  alen  messages  were  used  in  NEIC  alerts  in  1992.  In 
addition,  teleseismic  events  have  been  routinely  reevaluated  by  an  analyst,  and 
parameter  data  (onset  time,  amplitude,  period,  azimuth  and  slowness)  for  1 1649 
phases  (Le.,  32  phases  per  day)  have  been  sent  to  NEIC  (Golden,  CO)  and 
other  institutions  in  1992. 


In  1992,  the  data  acquisition  software  (TM  Science  Horizons)  was 
upgraded  several  times  and  became  robust.  The  RDAS-2(X)  appears  to  semi 
infrequently  data  with  corrupted  tixire  stamps.  These  data  appear  to  be  in  the 
future  (Hgure  1-4).  In  1992,  SHI  sireceeded  in  taking  care  of  this  problem  by 
software  inqnovements.  Table  1-3  shows  that  the  outage  due  to  future  data 
entering  the  data  acquisition  process  amounted  to  more  than  2  days  in  1992. 
One  old  problem  of  the  data  acquisitioo  strftware  has  still  remained  up  to  the 
end  c£  1992:  About  once  a  month,  the  SCHEDULER  driven  program  MRG 
crashes  2  consecutive  times.  This  results  in  a  data  loss  of  30  minutes.  Table  1-3 
shows  that  data  loss  due  to  this  im}blem  caused  outage  that  amounted  to  7J 
hours  in  1992. 


13  GERESS  Array  Hub 


The  GERESS  array  showed  stable  operetion  in  1992,  Table  1-2  shows  the 
xspfsxat  of  the  system  includittg  the  Bwhum  (X»ervaiory  (uptimes  directly 
conrespood  to  the  data  archived  in  Bochum).  Tire  avenge  uptime  was  96.9  %, 
ie.,  die  total  outage  amounts  to  113  days  in  1992.  The  hi^iest  uptime 
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Table  1*2 

GERESS  uptime  in  1992  (incL  Bochum  data  center) 


Month 

uptime 

% 

Jan 

99.7 

Feb 

99.0 

Mar 

87.8 

Apr 

91J 

May 

97.5 

Jun 

94.2 

Jul 

96.0 

Aug 

99.4 

Sep 

99.6 

Oct 

99.5 

Nov 

99.4 

Dec 

99.8 

Table  1*3 

Causes  of  GERESS  downtimes  in  1992  (incL  Bochum  data  center) 
(The  total  downtime  was  980300  sec) 


cause 

downtime 

downtime 

[sec] 

[%] 

Bochum  data  acq.  software  (SS2DL} 

16620 

1.7 

Bochum  data  acq.  software  (SCHEDULER) 

27000 

2.8 

Bochum  data  acq.  software  (future  data  received) 

178310 

18.2 

Bochum  data  acq.  software  (maintenance) 

4265 

0.4 

Bochum  data  acq.  woikstatioD 

97900 

10.0 

Bochum  QM 

2155 

02 

Bochum  power  outage 

19888 

2.0 

total  outage  Bochum 

346138 

353 

Communication  line 

9467 

1.0 

HUB  lAC  crashes 

329295 

33.6 

HUB  thunderstonn  damage 

203490 

20.8 

HUB  maintenance 

50719 

5.2 

HUB  power  ouci^ 

11226 

1.1 

total  outage  HUB 

594730 

60.7 

ocher 

29965 

3.0 
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(99.8  %)  was  observed  in  December  1992;  the  lowest  uptime  (87.8  %)  occinred 
in  March  1992  due'to  thundostorm  damage  at  the  array  sim.  Table  1-3  further 
specifies  the  causes  of  the  observed  outages.  The  causes  of  the  downtime  are 
clearly  related  to  the  well  known  problems  of  this  installation.  The  technical 
smtus  of  GERESS  is  summarused  in  chronological  order  in  Appendix  1-1 
(comprehensive  monthly  status  reports  are  available  upon  request). 

In  1992,  the  technical  installation  of  (jEKESS  could  not  be  considered 
fully  acceptable.  Problems  have  been  identified  during  the  GERESS  System 
Vexification  Test  in  September  1991  (Golden  et  aL,  1S191):  A  persistent  problem 
is  that  data  gaps  ate  inherent  in  the  lAC  The  nuntiier  of  I-second  gaps  per 
data  channel  were  found  to  form  4  groups.  These  4  groups  conespond  to  the  4 
communication  boards  (digiboatds)  on  the  array-controller  (Table  1-4).  The  rea¬ 
son  is  that  the  digiboards  perform  sporadic  self-resets,  loosing  about  9  seconds 
of  data  each  time  (Hgure  1-5).  Therefore,  the  overall  data  l<xss  is  significantly 
greater  than  the  loss  in  data  transmission  as  specified  and  observed  by  the  Ger¬ 
man  Telekom  (1  bit  in  10£-K)6). 

Based  on  the  findings  of  the  GERESS  System  Verification  Test,  Teledyne 
Geotech  started  corrective  actirni  under  warranty.  In  May  1992,  a  hardware 
(lAQ  and  software  upgrade  (LAC  and  RDAS-200)  was  sent  by  Teledyne  and 
installftri  by  Bochum  engineers.  The  idea  of  the  upgrade  was  to  make  Ae  baud 
rates  between  the  RDAS’es  and  digiboards  in  the  lAC  different.  Unfortunately, 
tins  upgrade  did  not  solve  the  data  gap  problem  (Table  1-4,  Figure  1-5).  At 
Dallas,  Teledyne  Geotech  has  continu^  to  work  on  this  problem;  but  at  the 
end  of  1992,  die  data  gap  probitsn  still  awaits  a  solution. 


Table  1-4 

Number  of  1  secmid  data  gaps  averaged  fix-  a  specific  digiboard  (DB).  Results 
of  6  tests  are  given  (times  were  adjusted  to  a  test  duration  of  241896  seomds). 


DB 

29.10.92 

21.8.92 

15.1Z91 

18.9.91 

5.6,91 

15.11.90 

1 

161 

137 

409 

522 

176 

274 

2 

270 

355 

511 

1084 

527 

459 

3 

131 

216 

181 

194 

240 

? 

4 

103 

158 

99 

161 

97 

? 

It  has  betel  noted  that  the  optical  data  links  between  RDAS'es  and  HUB 
are  very  sensitive  to  lightning  damage  due  to  destroctitm  of  optical  modems 
wh^  are  fed  by  DC  power  fiom  the  HUB.  This  means  that  the  optical  links 
have  to  undergo  mainrenanre  wmk  comparatively  often.  Since  proper  equip- 
ment  to  adjust  these  links  is  extremely  expensive  (and  does  n<»  exist  at 
GERESS),  die  adjustment  is  dcaie  by  trial  and  error.  A  deviatimi  or  drift  from 
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the  (^dmum  adjustment  means  that  time  synchronization  pulses  (1  pps  signals) 
can  be  distort^  hence  the  synchronization  of  the  array  is  fundazmntaily 
affected  An  example:  On  Nov  3,  it  was  noted  that  C5  was  out  of  synchroniza* 
don  by  0.4  seccnds.  The  exact  start  time  of  desynchronizadtm  was  deumnined 
as  Aug  17.  At  that  time,  strong  thunderstonm  had  damaged  various  channels  of 
the  array  repeatedly  and  also  damaged  the  opdcal  modem  of  CS.  Unfortunately 
this  had  not  been  noticed  at  that  dme.  One  problem  is  that  the  s^^hronizadon 
can  only  be  checked  by  looking  at  strong  teleseismic  onsets.  Another  indicator 
of  a  possible  misadjustment  in  the  optical  links  are  self-tesets  of  RDAS’es 
(e.g.;  D5  (Jan),  D8  (Apr),  B1  (May),  B3  (May),  B1  (Jun),  C7  (Jun),  C5  (Jun, 
Jul),  D8  (J;ai-Oct),  D9  (Aug),  D1  (Sep),  A2  (Sep-Dec),  B3  (Oct)).  Fot  stable 
channel  synchronizadra,  possibilities  to  provide  each  RDAS  with  an  indepen¬ 
dent  dme  receiver  should  be  invesdgat^  However  low  cost  GPS  receivers 
reach  a  synchronizadon  (+/-  1  msec)  that  seems  not  sufficient. 

The  cunent  release  of  RDAS-2(X)  software  installed  at  GERESS  shows  a 
problem:  A  stadon  does  not  synchronize  with  the  GERESS  master  clock  after  it 
comes  on-line  afnr  a  power  shutdown.  In  addidon,  the  RDAS  does  not  exe¬ 
cute  dxne-sets  piompdy,  sometimes  a  reinit  followed  by  a  dme-set  works, 
sometimes  ^veral  time-sets  are  necessary.  Furthomore,  various  channels  have 
shown  "artificial”  data  (i.e.,  pure  electronic  noise  that  mimics  seismic  noise 
pretty  closely  in  the  dme  domain  (e.g.  Figure  1-6)).  These  instances  appear 
after  a  power  up  or  after  a  self-reset  and  vanish  only  afttr  a  manual  reinit  fiom 
Bochum. 

Some  channels  of  GERESS  show  a  significant  higher  noise  level  due  to 
site  effects  (Jbst,  1993).  The  vaults  of  AO,  A3,  and  B3  could  not  be  placed  on 
bedrock  but  ate  <m  very  local  pockets  of  weathered  granite  or  gneiss.  In  order 
to  move  these  channels  within  the  ri^  restricdoos  set  by  the  array  layout,  a 
shallow  refiiacdon  study  was  perfonned  near  stadon  AO  and  A3  in  cooperadon 
mdi  die  BGR  (Gesiext^i,  1993).  This  was  the  first  step  of  a  mapping  project 
of  die  bedrock  around  the  cridcal  stadons.  Gestermaim  (1993)  conclud^  that 
bedrock  can  be  reached  at  a  depth  of  4S  m  below  AO  and  35  m  below  A3.  In 
the  vidnity  of  AO,  a  syncline  seems  to  be  present  and  a  lateral  shift  of  the  sta¬ 
don  location  is  pcohibidve.  Therefne,  a  bore-hole  seismometer  should  be 
emplaced.  For  A3  and  B3,  furthtx  mapping  of  the  basement  depth  should 
en^ie  a  lateral  shift  of  stadon  posidons. 


1.4  Condusioo 

At  die  GER..2SS  data  center  Bochum,  18  regional  evoits  have  been 
aimomarically  located  daily  on  average.  In  addidon,  extotsive  parameter  dam 
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for  32  teiesetsmic  phases  have  been  interacdveiy  ctetennined  each  day  on  aver¬ 
age.  The  GERESS  data  center  Bochum  is  the  only  German  staoOT  that  provides 
NEIC  with  alen  messages.  In  conclusion,  GERESS  is  the  roost  sensitive 
seismological  station  in  Central  Europe. 

The  Bochum  obsetvatoiy  routine  included  data  aichival  and  automatic  on¬ 
line  analysis.  The  automatic  locations  of  regional  events  render  a  good  first 
soludtm,  in  generaL  Since  October  1992,  regional  locations  ate  reevaluated  by 
analysts.  Results  of  the  interactive  re-analysis  of  d^a  have  ^fn 

used  by  NEIC  for  the  PDE.  In  addition,  interested  institutions  have  used 
parameter  data  of  GERESS  for  their  tasks.  The  Bochum  observatory  is  an 
qren  station  similar  to  the  German  Regional  Network  of  Broad-Band  Stations. 
In  addition,  Bochum  serves  as  data  center  where  various  bulletins  (e.g., 
automatic  NEIC  locations,  NEIC  alerts,  QED,  PDE  (weekly),  PDE  (monthly), 
fast  CMT’s  from  Harvard,  onset  list  of  YKA)  are  available  for  the  interested 
user. 

The  GERESS  data  center  Bochum  has  been  responsible  for  GERESS 
maintenance.  The  GERESS  uptime  which  includes  the  Bochum  observatory 
was  96.9  %  in  1992.  This  demonstrates  the  high  level  of  reliability  of  the 
installation  and  its  successful  maintenance.  On  the  other  hand,  the  technical 
installatkm  of  GERESS  cannot  be  considered  fully  acceptable  for  the  period 
covered  in  this  report  Problems  identified  during  the  GERESS  System 
Verification  Test  in  September  1991  (Golden  et  al.,  1991)  have  remained 
unsolved  in  1992. 
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Appendix  1-1:  Chronologicai  Status  List  GERESS 

Cafifanuion  sequence:  Jan-April  8:00-8:13;  May-Dec  03KX)-03;13.  D8  and  D9  show 
abnormal  amplitudes  at  calibration. 

Jan.  1992: 

-  Installation  of  test  version  of  the  Data  Request  Manager  (DRM)  by  Krake  Inc. 

Feb.  1992: 

-  lAC  crashes:  2 

-  Feb.  10  -  14,  maintenance  visit  by  RUB  engineer. 

-  Desynchronization  on  B1  (1  msec  due  to  defective  clock  board). 

-  RDAS-200  software  update  (version  3.01)  loaded  to  C2B  (desynchronization 
problem  BB13  on  a  RDAS-^IO  (Golden  et  aL.  1991)  was  not  fixed). 

•  Installation  of  hardware  for  power  cycling  the  lAC  remotely  (e.g.,  from 
Bochum)  (information  provided  on  the  PSAC  screen  of  each  RDAS-200  can  also 
be  remotely  obtained  for  diagnostics  at  Bochum). 

-  At  Bochum,  a  software  update  for  MONITOR_DL  (version  1.5J!,  SHI)  was 
installed  (several  problems  related  to  future  data  generated  on  the  data  acquisition 
system  remain). 

Mar.  1992: 

-  lAC  crashes:  2 

.  Data  txansmisskMi  stopped  foe  3  days  due  to  thunderstorm  damage  at  the  HUB. 

-  Mar  23  -  25,  emergency  maintenance  visit  by  RUB  engineen. 

•  RDAS-2(X)  software  update  (versiem  3.01)  loaded  to  C2B  (desynchronization 
problem  BB13  on  a  RDAS-200  (Golden  et  aL,  1991)  fixed). 

-  At  Bochum,  a  software  update  for  MONTTORJDL  (version  1.52,  SHI)  was 
installed  (several  problems  related  to  future  data  generated  on  the  data  acquisition 
system  remain). 

Apr.  1992: 

-  lAC  crashes:  5 

-  Apr  9-10,  maintenance  visit  by  RUB  engineer. 

-  Apr  23,  synchronous  self  reset  all  channels:  clock  problem. 

-  Officers  of  Bundesgrenzschutz  (BGS)  started  to  use  pan  of  HUB  as  cffice. 

May  1992: 

•  lAC  crashes:  8 

-  May  11-15,  maintenance  visit  by  RUB  engineers;  May  19-22  emergency  mainte¬ 
nance  visit 

-  change  of  band  rate  between  RDAS'es  and  lAC 

-  hardware  changes  lAC  (eeproms,  digiboards,  sdlc  board),  new  softwatt  lACexe 
(version  2.00) 

-  software  chrages  of  RDAS  (3.30) 
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•  digib(»id  2  failed. 

.  GERESS  DCF'77  clock  defective  and  tempraanly  repaned  by  RUB  enginears. 

•  vault  inspecdtM)  by  SAG  (sttbctmtractor  of  LahxDeyer  Inc.),  no  leaks  fcnuid. 

Jun.  1992: 

-  lAC  crashes:  2 

•  tfaunderstonn  damage  on  10  channels 

-  June  3-S  and  June  19*26,  maintenance  visits  by  RUB  engineers. 

-  zq>lacement  DCF-77  clock  installed  in  parallel  to  defective  one  interface  prob¬ 
lems). 

-  installation  of  new  digiboaid  2 

-  power  outage  on  HUB 

•  GERESS  GRAND  OPENING 

•  Warranty  repairs  M.  Browne  (Teledyne-Geotech) 

•  new  release  of  data  acquisiticm  softwm  (SHI):  DL2CB  (3.9.1),  CB2DB  (3.9), 
MONITOR_DL  (1.6). 

JuL  1992: 

.  LAC  crashes:  1 

.  thunderstorm  damage  cm  1  channel  (Jul  3),  15  channels  (Jul  10),  and  6  channels 
(Jul22) 

•  July  13-14  and  July  27,  emergency  and  regular  maintenance  visits  by  RUB 
engineers. 

•  new  release  of  data  acquismop  software  (SHI):  SS2DL  (2.13), 
MONITOR_PRGS  (13.1),  MRG  (2.10),  RUNMRG.  Severe  malfunction  of  MRG 
(2.10)  documented. 

-  internet  number  for  data  analysis  workstation 

-  upgrade  of  data  request  managa  (DRM)  by  Krake  Inc. 

Aug.  1992: 

-  lAC  crashes:  1 

-  thunderstorm  damage  of  2  channels  (Aug  10),  11  channels  (Aug  21),  and  2 
channels  (Aug  26);  field  nip  not  necessaiy. 

-  extensive  data  gap  testing  under  tbe  direction  of  Teledyne. 

-  new  release  of  dm  acquisition  software  (SHI):  MRG  (2.10.1),  RUNMRG. 

-  power  outage  HUB 

-  power  glitch  at  Bochum  data  center  due  to  dinnderstorm 

Sep.  1992: 

-  lAC  crashes:  2 

•  tfaunderstonn  damage  on  1  channel  (field  trip  not  necessary). 

-  data  quality  received  at  data  center  in  Norway  itetetiatated  to  the  pmnt  duu  data 
are  useless  (Sep  21).  This  conelates  in  time  to  extended  outages  on  digtboard  2 
The  GERESS  data  received  at  data  center  in  Kjeller  remained  unusable  to  the  end 
of  1992. 
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Oct  1992: 

-  diunderstonn  damage  of  10  sianons  (Oct  7). 

•  Oa  5>24,  mainttnance  ivock  at  HUB. 

•  installation  of  in^noved  thundexsmrm  protection  equipment  and  power  condi- 
dmiing  unit  at  each  RDAS 

-  loading  new  RDAS  configuradon  files  (calibration  sequence  included) 

•  invesdgating  calibration  protdems  on  D7  irtpuxed),  D8,  D9. 

•  installadon  of  new  DCF-TJ  aasar  clock  for 

•  Connecdmi  of  Saeckeisen  (STS-2)  at  C2  (80  Hz)  on  Oct  8  (disconnection  of 
BB>13  and  GS-I3  @  120  Hz). 

>  start  of  reviewing  regitmai  event  locations  by  analysts 

Nov.  1992: 

Nov  25'26,  maintenance  work  at-HUB. 

>  supervision  of  digging  near  nutin  cable  route  to  HUB 

Dec.  1992: 

Dec  14-18,  maintenance  work  at  HUB. 

-  mstallarion  of  improved  dmnderstorm  protection  equipment  at  HUB. 

•  BB-13  (vertical)  reconnected 

•  absolute  n™  displayed  on  DCF-77  master  clock  and  in  the  lAC  different  by 
minutes 

•  Bochum:  upgrade  (versicm  1.19)  data  request  manager  (DRM  by  Kxake  lix;.) 
installed 
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Hgure  1-1:  GBRESS  Observatory  Bochum.  Data  Flow. 
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GERESS  BOCHUM:  AUTOMATIC  LOCATIONS  FOR  1992 


Figure  1-2 
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GERESS  BOCHUM:  ANALYST  REVIEWED  i  nrATIDNS  nrT.npr  i 


Figure  1-3 
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Figure  1*5:  QERESS  wovefomi  displuy  for  nil  chnnnels  on  December  II,  |992  11:42:00. 
Noie  liie  "9-second''  dnin  gaps  on  all  channels  of  digiboard  3  (B3  is  out). 


2.  REPORT  ON  GERESS  INSTRUMENT 

CALIBRATION 


Jan  Wuster 


2.1  Abstract 

GERESS  local  magnitudes  have  been  observed  to  be  systematically  low  by  at 
least  0.5  magnitude  units  in  comparison  with  neighbouring  stations.  In  order 
to  exclude  the  possibility  that  this  effect  is  due  to  some  sort  of  instrumental 
or  procedural  error,  a  thorough  check  on  GERESS  instrument  calibration 
was  performed  in  October  and  November  1992.  Components  of  the  sys¬ 
tem  which  were  checked  include:  GS13  generator  constant,  transfer  function 
and  calibration  coil  motor  constant,  RDAS  calibration  signals,  RDAS  pre- 
ampliffcation  and  A/D-conversion.  Although  several  minor  deviations  from 
technical  specifications  were  found,  none  of  them  can  atccount  for  a  system¬ 
atic  diminution  of  amplitudes  by  a  factor  of  2  or  more.  Consequently  the  low 
magnitudes  need  to  be  explained  otherwise. 

2.2  The  Problem 

At  GERESS,  Richter  local  magnitudes  are  determined  from  Wood- Anderson 
simulated  traces  (seee.g.  Jost,  1992)  using  an  amplitude-distance  relationship 
(or  Western  Europe  (Nicolas,  1982.)  To  take  advantage  of  the  improved 
signal/noise  ratio,  standard  values  reported  in  the  GERESS  event,  lists 
are  calculated  from  the  vertical  beam  trace.  Horizontal  beams  and  a  single 
horizontal  station  (GEC2)  have  also  been  used  for  comparison.  The  signal 
loss  due  to  beam-forming  was  found  to  be  0.2  magnitude  units  on  average 
(Schweitzer,  pers.  comm.)  Even  if  this  is  taken  into  account,  GERESS  Aft- 
values  are  systematically  low  in  comparison  with  neighbouring  stations  and 
netwwks. 

The  effect  has  been  observed  for  some  time  by  GERESS  and  GRF 
staff,  it  was  studied  quantitatively  by^  Kradolfer  (1992),  who  found  that  the 
Aft-values  from  GERESS  were  on  the  average  0.8  ±  0.3  magnitude  units 
lower  than  the  corresponding  values  calculated  by  the  Swiss  Seismoiogical 
Service  for  a  sample  of  42  Swiss  events  from  Nov  19  1991  to  May  18  1992, 
although  part  of  tMs  difference  could  be  attributed  to  the  different  procedures 
used  for  local  magnitude  determination,  ms-determinations  at  GERESS  for 
tdieseismic  events  also  tend  to  be  lower  than  the  values  reported  by  NEIC 
(Schweitzer,  this  report,  chapter  3) 


18 


Before  attributing  these  findings  to  external  factors,  we  must  make 
sure  that  any  internal  causes  can  be  excluded.  Therefore  it  was  decided  to 
check  thoroughly  all  relevant  components  of  the  GERESS  signal  transmission 
system  for  compliance  with  technical  specifications. 


2.3  Routine  calibration  checks 

Routine  calibration  checks  are  being  performed  every  day  between  03:00  and 
03:13  UT.  A  calibration  signal  (4500  mV  zero-to-peak  nominal  amplitude, 
1  Hz  frequency,  50  s  duration)  is  being  fed  into  all  channels  following  a  sched¬ 
ule  published  most  recently  in  the  GERESS  status  report  of  August  ’92  (Jost 
1992.)  [Note  that  4500  mV  is  the  nominal  voltage  value  downloaded  to  the 
RDASs.  For  actual  voltage  values  see  below  under  2..5.]  In  each  seismometer 
the  calibration  signal  drives  a  calibration  coil  in  series  with  a  capctcitor,  mak¬ 
ing  the  electrically  induced  velocity  of  the  seismometer  mass  prop>orti6nal  to 
the  applied  voltage  (Teledyne  1988,  page  628/19.)  Teledyne  measured  the 
transfer  fimction  ^(s)  of  the  electrically  driven  seismometer  (Horn  1990) 
and  express  them  in  terms  of  poles,  zeros  and  gain  [the  gain  being  valid  for 
the  nominal  calibration  voltage  values.]  In  theory,  transfer  functions  of  all 
GERESS  GS  13-channels  should  be  identical,  because  effective  generator  con¬ 
stant,  natmal  frequency  and  effective  damping  have  been  normalized  to  2000 
Vs/m,  1  Hz  and  0.75  critical  damping  respectively.  The  amplitude  of  the 
transfer  function  ev^uated  at  I  Hz  is  H{2Ti  •  IHz)  «  70727.3  ctmnisfV oH. 
This  means  that  with  Ueai  —  4.5  V  we  expect  to  measure  318273  digital 
counts  zero-to-peak.  Table  2  -1  shows. the  results  of  a  comparison  of  ampli¬ 
tude  values  measured  during  a  routine  calibration  on  Oct  28  1992  with  the 
theoretical  values,  giving  the  deviation  in  percent. 

13  of  33  short  period  channels  are  within  specifications  (±1%),  another 
14  channeb  are  within  ±5%  of  the  theoretical  value.  Deviation  is  not  ac¬ 
ceptable  for  channeb  D4jz  (—5.1%),  Dl^  (—6.4%),  D7j5z  (—7.8%)  and 
Dl.5n  (-t-12.2%).  Values  for  DSjsz  (factor  26.4)  and  D9.sz  (factor  1/185.4) 
are  drastically  off.  These  drastic  deviations  axe  not  observed  in  seismic  sig¬ 
nal  or  noise  amplitudes  and  will  be  considered  sq>arately  below  in  section 
2..10.  The  percentage  values  in  Table  2!-l  were  checked  on  several  days  and 
were  found  to  be  quite  stable.  Daily  fluctuations  t}rpically  affect  the  second 
decimal  place  only. 

It  is  concluded  thid  although  some  channeb  axe  definitely  outside  spec¬ 
ifications,  GS13  overall  channel  sensitivity  to  electrical  excitation  at  1  Hz  is 
as  expected. 
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Table  2-1:  Check  on  Calibration  of  GS13  40  Hz  channels  at 
1992-302:03.00.00.000 

Vertical  Channels,  CALNORM  •  318273 


1 

GEAO.sz  316598  « 

-0.526  X 

2 

GEAl.sz  320279  - 

0.630  X 

3 

GEA2.8Z  311520  « 

-2.12  X 

4 

GEA3.SZ  320545  - 

0.714  X 

5 

GEBl.sz  319350  - 

0.338  X 

6 

GEB2_sz  316753  - 

-0.478  X 

7 

GEB3.SZ  309209  « 

-2.85  X 

8 

GEB4.SZ  319951  « 

0.527  X 

9 

GEB5.SZ  310433  > 

-2.46  X 

10 

GECl.sz  315426  - 

-0.895  X 

11 

GEC2.SZ  312798  » 

-1.72  X 

12  GEC3.SZ  315587  » 

-0.844  X 

13 

GEC4.SZ  311976  » 

-1.98  X 

14  GEC5.SZ  311985  « 

-1.98  X 

15 

GEC6.SZ  316194  > 

-0,653  X 

16  GEC7.SZ  304366  • 

-4.37  X 

17  GEDl.sz  297857  « 

-6.41  X 

18  GED2.8Z  325289  • 

2.20  X 

19  GED3.SZ  317074  « 

-0.377  X 

20 

GED4.SZ  302023  « 

-5.11  X 

21 

GED5.8Z  307657  • 

-3.34  X 

22  GED6.SZ  318786  • 

0.161  X 

23  GED7.8Z  293462  « 

-7.80  X 

24  GEDS.sz  0.838861E+07  — —  too  big  for  cal 

25  GED9.8Z  1717  —  too  snail  for  cal 

Horizontal  Channels,  CALNORM  ■  318273 

26  GEA2.sn  318578  -  0.096  X 

27  GEA2.se  315817  •  -0.772  X 

28  GEDl.sn  357015  -  12.17  X 

29  GEDl.se  314265  -  -1.26  X 

30  GED4.sn  312781  •  -1.73  X 

31  GED4_se  312670  «  -1.76  X 

32  GED7.sn  309569  -  -2.73  X 

33  GED7.se  308426  •  -3.09  X 
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Figure  2-1;  Sinusoid  signals  of  various  frequencies  and  fixed  axnplitud*^  were 
applied  to  the  calibration  coil  on  channel  GED2  jz: 


2.4  'Bransfer  functions 

The  GS13  sensitivity  for  ground  velocity  has  a  plateau  for  frequencies  >  1 
and  a  descending  flank  for  frequencies  <  1  Hz,  the  natural  frequency  of 
1  Hz  being  the  3  dB  point.  The  exact  shape  of  the  sensitivity  curve  in  the 
transition  region  is  relevant  for  local  magnitude  determination,  because  the 
largest  amplitudes  in  the  S-wavetrain  often  have  periods  around  T  =  I  s. 

An  extended  electrical  calibration  was  performed  on  two  randomly  se- 
,  lected  vertical  channeb  (D2_sz  and  B5jz.)  To  this  end  sinusoidal  signals  of 
various  frequencies  in  the  range  0.1 ...  I  Hz  and  a  constant  amplitude  were 
j4>plied  to  the  calibration  coils  of  the  respective  seismometers  (Fig  2-1.) 
The  signal  amplitudes  measured  match  the  amplitude  of  the  transfer  func¬ 
tion  quite  well  (see  Fig  2-2),  except  in  the  case  of  the  calibration  signal  of  0.1 
Hz,  which  is  obviously  not  produced  correctly  by  the  RDAS  signal  generat¬ 
ing  unit  (Figure  2-3.)  The  small  sensitivity  surplus  at  higher  frequencies  on 
Channel  D2^  is  best  brought  out  by  a  half-logarithmic  plot  (Figure  2-4.) 
Although  measured  values  lie  on  the  theoretical  curve  (the  solid  line) 
at  and  below  1  Hz,  they  are  significantly  higher  above  1  Hz.  To  adjust 
the  theoretical  curve  to  the  measured  value  at  5  Hz,  the  static  gain  would 
needs  raised  by  6%  (the  dashed  line),  at  the  expense  of  a  poor  fit  in  the  1 
Hz-regi<m.  The  second  channel  checked  was  found  to  be  within  specification 
over  the  whtde  frequency  range  covered  (Figure  2-5.) 
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Figure  2-4:  D2jsz  calibration  measurements  and  theoretical  transfer  function 
amplitude:  half-logarithmic  scale 


xio» 


It  is  concluded  that  even  when  sensitivity  matches  the  specified  value 
closely  at  1  Hz,  this  may  not  necessarily  be  true  at  other  frequencies.  But 
the  1  Hz-re^on  is  most  important  for  Mi  determination  and  hence  the 
observed  amplitude  deviations  at  higher  frequencies  cannot  be  responsible 
for  a  systematic  diminution  of  magnitudes. 

2.5  RDAS  signal  generator 

So  far,  channel  sensitivity  has  been  checked  using  electrical  excitation  by 
the  RDAS  signal  generating  unit.  But  this  unit  itself  might  be  faulty  or 
the  electrical  excitation  might  for  some  reason  differ  from  seismic  excitation. 
Therefore  an  effort  was  made  to  test  all  the  components  which  determine  the 
amplitude  sep2irately,  starting  with  the  signal  generator. 

Here,  indeed,  a  major  deviation  was  found:  Calibration  signals  output 
by  the  RDAS  were  displayed  on  an  oscilloscope  in  vaults  D7,  D8  and  D9 
at  various  voltage  settings  from  10  mV  to  4500  mV  and  were  consistently 
found  to  be  too  high  by  a  factor  of  2  (see  Table  2!-2.)  If,  for  example, 
the  calibration  voltage  has  been  downloaded  to  be  100  mVop,  the  actual 
signal  output  will  have  an  amplitude  of  400  mV^,  equivalent  to  200  mVop. 
In  order  to  explain  the  fact  that  overall  sensitivity  to  electrical  excitation 
has  been  found  to  be  in  conformity  with  the  calibration  transfer  function 
supplied,  while  calibration  signal  amplitudes  are  obviously  in  error,  it  must 
be  assumed,  that  the  calibration  transfer  function  has  been  originally  derived 
using  the  excessive  RDAS  signal  amplitudes.  Electrical  excitation  without 
use  of  the  RDAS  signal  generator  is  considered  below  in  section  2.. 9. 

A  side  result  of  the  checks  on  the  signal  generators  is,  that  the  pecular- 
Hies  observed  in  routine  calibration  (see  section  2..3)  on  channels  D8-sz  and 
D9jsz  are  not  caused  by  the  signal  generators. 


2.6  Pre-amplifier  and  A/D-converter 

Pre-amplification  and  A/D-converting  stages  were  not  checked  separately  but 
as  a  unit.  Calibration  signals  of  known  amplitudes  were  fed  directly  into  the 
RDAS  signal  input  ports  on  channels  D8.az  and  D9..SZ,  the  seismometer  being 
disconnected.  Amplitudes  measured  on  the  recorded  trace  (in  digital  counts) 
were  converted  to  signal  input  levels  using  the  specified  A/D-converter  LSB  =: 
2.3841857  itVj count  and  the  pre-amplification  factor  of  V  s  31.6228  (30  dB.) 
The  results  are  shown  in  Table  2!-2.  Calculated  input  signal  leveb  (zero- 
to-peak)  are  close  to  the  signal  amplitudes  measured  with  an  oscilloscope 
(zoo-to-peak),  as  expected.  Oscilloscope  measurements  are  only  exact  up  to 
dbl0%,  but  if  we  assume,  that  the  output  of  the  signal  generators  is  exactly 
twice  the  value  specified  (see  above  under  2!-2)  there  is  less  than  2%  deviation 
between  nominal  amplitudes  times  2  and  measured  amplitudes  (assuming 


24 


Table  2-2:  Calibration  signal  amplitudes  at  1  Hz 


channel 

nominal  amp 
(mVbp) 

measured  amp 
[mVop]  ±  10% 

max  trace  amp 
dig.  counts  op 

calc,  input  amp 
[mVop] 

D9-SZ 

25 

2.612  •  10^ 

19.693 

40 

5.262  •  10® 

39.673 

85 

1.048  •  10® 

79.014 

60 

120 

1.571  •  10® 

118.45 

80 

165 

2.097 . 10® 

158.10 

100 

210 

2.651  •  10® 

199.87 

D8jsz 

80 

79.104 

120 

118.603 

160 

158.200 

100 

200 

2.6523  •  10® 

199.969 

specified  LSB  and  V.)  The  fact  that  all  deviations  are  negative  indicates, 
that  either  the  signal  generator  output  level  multiplied  by  2  or  LSB  and  V 
are  about  1%  lower  than  expected  on  the  two  channeb  checked. 

It  is  concluded  that  the  pre-amplification  and  A/D-conversion  stages 
of  the  RDAS  are  not  responsible  for  a  diminution  of  magnitudes.  Also,  cali¬ 
bration  irregularities  observed  on  chaimels  D8-sz  and  D9-sz  (section  2.. 3)  are 
not  due  to  these  stages,  they  were  reproduced  after  changing  the  respective 
RBAS  A/D-converter  card. 


2.7  Seismometer  generator  constant 

Weight-lifting  tests  were  performed  on  channels  D7..SZ,  D8jz  and  D9^z.  The 
results  could  only  be  used  qualitatively  to  confirm  that  generator  constants  at 
D8  and  D9  (showing  irregularities  in  calibration)  do  not  drastically  differ  from 
D7  (not  showing  irregularities.)  An  indirect  method  was  used  to  confirm  the 
specified  GS13  generator  constant  of  2000  Vs/m.  An  STS2-seismometer  with 
a  specified  generator  constant  of  1500  Vs/m  has  been  temporarily  installed 
in  the  GEC2B  vault.  Its  3-component  output  is  recorded  on  the  C2Jm,  Jiz 
and  Jie  high  frequency  channels,  and  the  STS2-C2_hz  trace  can  be  compared 
to  the  GS13-C2.SZ  trace  measured  at  the  GEC2A  vault.  The  STS2  transfer 
function,  however,  is  flat  to  velocity  up  to  a  period  of  T  s  120s,  so  the  traces 
cannot  be  compared  directly  (Figure  2-6.)  But  from  each  of  the  different 
traces  a  Wood/Anderson  trace  can  be  simulated  using  suitable  HR  filters 
constructed  on  the  basis  of  the  specified  transfer  functions  and  generator 
constants. 

Figure  2-6  shows  the  maximum  S-wavegroup  amplitudes  of  an  event 
from  Switzerland,  from  which  Ml  would  be  determined.  The  upper  2  traces 
are  the  original  STS2  and  GS13  traces,  the  lower  2  traces  are  the  simulated 
Wood/ Anderson  traces.  The  latter  are  indistinguishable  in  signal  shape,  but 
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Figure  2-6:  Time-mndow  containing  the  S-wavegroup  maximum  amplitude 
of  an  event  from  Switzerland.  Upper  two  traces:  STS2  and  GSI3  origi¬ 
nal  seismograms.  Lower  two  traces:  Wood/Anderson  traces  simulated  from 
STS2  and  GS13,  respectively:  the  lower  traces  are  indistingmshabie  except 
for  a  7%  differoice  in  static  gain. 
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Figure  2-7:  GS13  sensitivity  measured  on  shaking  table  without  output  net¬ 
work 


there  is  a  diffoence  in  measured  zero-to-peak  maximum  amplitude  in  mm, 
the  trace  simulated  from  STS2  ]delding  a  value  approximately  7%  higher. 
Even  before  the  substitution  of  the  STS2  for  the  original  GS13  instrument 
on  C2Jiz,  this  channel  yielded  approximately  10%  higher  amplitude  values 
in  seismograms,  an  effect  which  was  not  observed  in  electrical  calibration.  It 
is  concluded  that  the  performances  of  STS2  and  GS13  reciprocally  support 
each  other  and  that  the  generator  constants  of  the  two  instruments  are  within 
10%  of  specifications. 


2.8  Shaking  table  experiments 

The  vertical  GS13  seismometer  with  serial  No.  181,  formerly  deployed  in  vault 
GEA2,  was  shipped  to  Bochiim  and  subjected  to  some  shaking  table  tests. 
Setup  and  performance  of  these  experiments  followed  Hinzen  (1979).  By 
means  of  the  shaking  table  it  is  possible  to  determine  seismonaeter  output  as  a 
function  of  gnmnd  movement,  i.e.  perform  a  seismic  calibration.  In  our  case, 
output  voltage  maximum  amplitude  was  set  in  relation  to  ground  velocity 
maximum  amplitude.  This  velocity  sensitivity  has  the  unit  V/{m/s)  &=  Vs/m 
and  is  identical  to  the  seismometer  generator  constant  in  the  plateau  well 
above  the  natural  frequency. 

In  the  first  series  of  tests  the  sosmometer  coil  output  was  measured 
directly,  not  using  the  resistor  network  normally  employed  (Figure  2-7.)  The 
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Figure  2-8:  GS13  sensitivity  measured  on  shaking  table  with  output  network 


system  is  underdamped,  but  values  to  the  right  of  1  Hz  lie  well  on  the 
theoretical  curve  calculated  using  the  GS13  transfer  function  and  the  '^native" 
generator  constant  of  the  instrument  No.  181,  Gk  ~  2316  Vsfm  (Horn.  1990). 
Note  that  the  broad  resonance  at  57  is  not  observed. 

The  function  of  the  output  resistor  network  is  to  normalize  both  effec* 
tive  damping  K  and  effective  generator  constant  Gt  to  convenient  values  for 
all  instruments  in  a  set  (Teledyne  1988,  page  628/19.)  According  to  Horn 
(1990),  the  effective  values  for  all  GER^S  GS13  seismometers  are  set  to 
Xr  s  0.775  '  critical  and  G«  =  2000  Vsfm  respectively.  In  the  second  series 
of  tests  the  seismometer  coil  was  connected  via  the  output  resistor  network. 
Figure  2>8  shows  that  the  measiued  values  lie  significantly  below  the  theo¬ 
retical  curve  calculated  using  G  =  2000  Vajm,  A  reasonable  fit  is  achieved 
assuming  G  b  1800  Vs/m  as  shown  in  Figure  2-9  in  linear  scale.  Neither 
the  substantial  resonance  at  57  Hz  nor  the  sharp  drop  at  80  Hz  are  observed. 
There  are  hints  of  small  resonances  at  S.and  70  Hz^  but  this  question  should 
be  addressed  in  more  detail  after  some  imperfections  of  the  experimental 
setup  have  been  removed. 

It  is  concluded,  that  the  sensitivity  to  ground  velocity  of  the  seismome¬ 
ter  under  study  follows  the  amplitude  of  the  transfer  function  specified  by 
the  manufacturer  below  20  Hz.  Above  20  Hz  the  senritivity  remains  fiat 
to  vdodty  contrary  to  the  transfer  function  supplied  by  the  manufacturer, 
whidi  has  a  broad  resonance  p«k  between  30  and  70  Hz.  The  measured 
effective  generator  constant  is  10  %  below  spedfication. 
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Figure  2-9:  GS13  transfer  function  amplitude  fitted  to  measured  sensitivity 
values 


2.9  Electrical  excitation 


In  order  to  prove  the  hypothesis  on  the  calibration  transfer  functions  in  sec¬ 
tion  2..5,  the  functionality  of  the  calibration  coil  was  tested  independently  of 
the  RDAS  calibration  signal  generator.  The  calibration  coil  is  driven  in  series 
with  a  2  nF  capacitor,  which  makes  the  velocity  of  the  seismometer  mass  in¬ 
duced  by  the  calibration  coil  proportional  to  the  applied  calibration  voltage. 
The  proportionality  constant  calculated  according  to  Teledyne  (1988,  page 
628/10)  is 


C  •  M 

VELSENS  = 

m 


2.0  •  4.5  f  ^  ^  nm/s 


with  the  capacity  C,  the  calibration  coil  motor  constant  Me  and  the  sus¬ 
pended  mass  m. 

Now  the  seismometer  can  be  excited  applying  sinusoidal  signals  of  var¬ 
ious  frequencies  with  an  amplitude  of  Etti,  seismometer  output  amplitudes 

are  measured.  Ratios  V^IUed  thus  determined  can  be  compared  to 
thewetical  ratios  obtained  using  constant  VELSENS  and  the  seismometer 
transfer  frmction  on  one  hand  and  the  calibration  transfor  function  on  the 
other  hand.  Figure  2-10  shows,  that  the  measured  values  follow  the  first 
curve  (dadied  line)  reasonably  well,  while  the  latter  curve  (solid  line)  is 
definitively  too  hi{^  by  a  factor  of  2.  Note  that  again  the  resonance  at  57 
Hx  is  not  obsCTved. 
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Figure  2-10:  GS13  electrical  excitation,  voltage  ratio  measured  and  calcu* 
lated 


Thus  it  has  been  confirmed  that  the  calibration  transfer  function  sup¬ 
plied  by  the  manufacturer  must  only  be  used  in  connection  with  the  excessive 
RDAS  cal  signal  amplitudes.  The  problem  must  be  bom  in  mind  in  case  the 
RDAS  is  used  to  calibrate  a  difierent  seismometer,  or  the  GS13  is  calibrated 
without  the  RDAS.  This  peculiarity  has  no  effect  on  earthquake  magnitude 
determinations. 

2.10  Calibration  Irregularities 

The  channels  DSjsz  and  D9j5Z  were  found  to  be  drastically  off  in  routine  cal¬ 
ibrations  (see  above  2..3.)  RDAS  signal  generators,  pre-amplifiers  and  A/D- 
converters  could  be  excluded  as  sources  of  the  problem  (above  under  2.. 5  and 
2..6.)  In  the  case  of  D9  the  problem  (amplitudes  too  low)  disappeared  when 
the  seismometer  was  exchanged  with  the  one  previously  deployed  at  A2^. 
When  opportunity  arises,  the  sdsmometer  should  be  checked  at  Bochum  and 
sent  to  the  manufacturer  for  repair,  if  no  fix  can  be  found.  In  the  case  of 
D8  (amplitudes  too  high)  the  problem  persisted  even  after  the  exchange  of 
the  seismometer.  Here  the  reason  could  be,  that  the  RDAS  does  not  make 
the  necessary  internal  connection,  which  corresptmds  to  a  jumper  connec¬ 
tion  between  pins  H  and  U  of  the  GS13  connector  required  in  calibration 
mode  (Tetedyne  1988,  page  828/10.)  In  this  case  the  c^bration  coil  would 
be  driven  directly  by  the  signal  generator  without  a  capacitor  in  series.  An 
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exchange  of  the  appropriate  RDAS  board  should  fix  the  problem.  The  ir¬ 
regularities  described  occur  in  calibration  mode  only  and  are  irrelevant  for 
magnitude  determination. 


2.11  Conclusions 

The  following  is  a  resume  of  findings  encountered  in  the  check  on  GERESS 
instrument  calibration: 

1.  Only  13  of  33  channels  are  within  1%  of  specification  on  overall  sen¬ 
sitivity  at  1  Hz,  another  14  channels  are  vrithin  5%.  4  channels  are 
seriously  off,  2  channels  are  drastically  off.  In  the  case  of  the  latter, 
the  problem  is  restricted  to  calibration  mode. 

2.  The  RDAS  calibration  signal  generator  produces  sinusoidal  signals  with 
twice  the  nominal  amplitude.  Also,  it  cannot  produce  frequencies  below 
0.2  Hz. 

3.  Even  though  the  calibration  at  1  Hz  is  within  specifications,  a  deviation 
of  4-6%  has  been  observed  at  higher  frequencies  in  one  out  of  2  channels 
checked. 

4.  In  two  spot-checks,  LSB  and  pre-amplification  factor  were  found  to  be 
within  2%  of  specification. 

5.  The  effective  generator-constant  of  the  GS13-seismometer  at  GEC2.SZ 
was  found  to  be  within  10%  of  specification  by  comparison  with  an 
STS2  seismometer,  which  was  operated  in  parallel. 

6.  In  a  shaking- table  experiment  with  the  GSlS-seismometer  formerly  in¬ 
stalled  at  GEA2jsz,  the  effective  generator  constant  Gt  was  found  to 
be  10%  below  specification,  although  the  "native”  generator  constant 
Gk  was  foxmd  to  be  as  specified. 

7.  In  an  experiment  driving  the  seismometer  with  a  voltage  applied  to  the 
calibration  coil  in  series  with  a  capacitor,  the  sensitivity  of  induced  ve¬ 
locity  to  voltage  applied  was  found  to  be  as  specified.  The  calibration 
transfer  function  supplied  by  the  manufacturer  gives  excessive  sensi¬ 
tivities  by  a  factor  of  2.  This  effect  counteracts  the  elevated  RDAS 
calibration  signals  mentioned  above. 

8.  The  resonance  between  50  and  60  Hz  and  the  sharp  drop  above  70  Hz 
in  the  amplitude  of  the  GSl3  transfer  function  given  by  Teledyne  are 
not  observed.  Instead,  the  plateau  extends  to  100  Hz. 

None  of  these  findini^  justifies  the  assumption,  that  amplitudes  are 
i^tematically  diminished  by  the  GERESS  signal  processing  system  to  a 
degree  necessary  to  explain  the  low  local  magnitudes. 


31 


Acknowledgements 


Field  measurements  at  the  vaults  were  performed  by  Bernhard  Klotz  and 
Lothar  Kuhne.  The  latter  and  Andreas  Muller  assisted  with  the  shaking 
table  experiments.  Their  help  is  gratefully  acknowledged. 

References 

Hinzen,  Klaus-Gunter  (1979) 

Untersuehungen  zum  Antoendungsbcnich  einer  Schuttdtischanlage  bei 
der  Kalibrierung  von  Seismometem  und  Geophonen  unter  Verwendung 
verschiedener  Anregungsarten.  Diplomarbeit,  Institute  of  Geophysics, 
Ruhr-University  Bochum,  Germany,  July  1979. 

Horn,  Charles  (1990) 

CD  station  channel  transfer  characteristics.  Undated  and  unsigned 
manuscript  contained  in  a  letter  by  Charles  Horn  (Teledyne  Geotech)  to 
H.’P.  Harjes  (Institute  of  Geophysics,  Ruhr-University  Bochum)  dated 
Sept  11  1990.  The  letter  also  contains  data  sheets  on  individual  GER- 
ESS  instruments. 

Jost,  Michael  L.  (1992) 

Current  status  and  results  of  the  GERESS  daia  center  in  Bochum.  Pro¬ 
ceedings  of  the  GERESS  symposium,  Waldkirchen,  June  22-24  1992. 

Kradolfer,  Urs  (1992) 

Recent  seismicity  in  Switzerland  as  located  by  the  Swiss  network  and  the 
GERESS  array.  Proceedings  of  the  GERESS  symposium,  Waldkirchen, 
June  22-24  1992. 

Nicolas,  M.,  Massinon,  B.,  Mechler,  P.  &  Bouchon,  M.  (1982) 

Attenuation  of  regional  phases  in  Western  Europe.  Bull.  Seis.  Soc.  Am. 
72  2089-2106. 

Teledyne  Geotech  (1988) 

Operation  and  Maintenance  Manual.  Portable  Short-Period  Seismome¬ 
ter,  Model  GS-IS.  Stock  No.  990-55400-9800,  Teledyne  Geotech,  3401 
Shiloh  Rd.,  Garland,  Texas  75041,  USA. 


32 


3.  TELESEISMIC  DETECTION  AND  LOCATION 
CAPABILITIES  OF  THE  GERESS  ARRAY 


Johannes  Schweitzer 


3.1  Introduction 

Since  evaluation  of  the  GSE1T-2  experiment  (Schweitzer,  1992; 
Sdiweitzer  et  aL  1992)  in  1991,  it  is  known,  that  GERESS  is  the  most  sensitive 
seismological  station  in  Central  Europe.  In  October  1991,  the  GERESS  group  in 
Bodhum  started  to  reanalyse  teleseismic  ons^  d^ected  automatically  with  the 
GERESS  array.  The  reanalysed  parametm  were  sent  to  international  data  centers 
(NEIC,  EMSQ  and  to  any  interested  seismic  station.  The  reanalysis  was  neces¬ 
sary  because  the  automatically  picked  onset  times,  phase  identifications,  and  the 
res^  of  the  automatic  fk-analysis  (ray  parameters  and  azimuth's)  were  often 
errmeous. 

Since  October  18,  1991,  the  GERESS  onsets  were  associated  to  evoits  by 
NEIC  as  seen  in  die  weekly  or  monthly  Earthquake  Data  Report's  (EDR's)  of 
NEIC.  An  evaluation  of  GERESS  corttribution  gives  some  informations  about  the 
teleseismic  detection  curability  of  the  array  and  the  bias  of  observed  parameters 
sudi  as  aamuth,  ray  parameter,  onset  time,  and  amplitude. 

hi  die  following,  the  GERESS  contribution  to  the  EDR’s  (monthly  listings) 
from  October  18, 1991  to  August  1992  arul  to  the  EDR's  (weekly  listings)  from 
Sqstendrar  1992  to  Fdmiary  1993  is  evaluated.  Fig.  3-'  shows  all  events  for 
adikh  NEIC  associated  a  GERESS  {diase.  The  CARESS  oxitributions  can  be 
compared  widi  those  of  other  arrays  or  stations  adikfa  ooittiibated  routinely  to 
NEIC.  In  this  study,  the  contributions  of  the  Yellowknife  array  (YKA)  in  North¬ 
ern  Canada  and  die  contribution  of  the  singfe  seismic  station  Kaipmki^  Hory 
(KHC)  located  in  the  Bohemian  Massif  in  the  Czech  Republic  in  the  vicinity  (SS 
bn)  of  GERESS  were  used  comparison. 


3JI  GERESS  detection  capabfltties 

Every  event  in  the  range  between  and  IS(y*  epicentral  distanoe  ftv 
whidi  NEIC  reported  a  body  vi^e  magnitude  m;,  was  defined  as  an  obsovable 
tdeaeismic  event  (Note  drat  for  each  station  under  considenoioa,  the  observiMe 
events  are  different)  On  average,  CSERESS  reported  44.6  %,  YKA  59 J  %,  and 
KIK  26 J  %  of  dieae  events.  The  differences  between  GERESS  and  KHC  can  be 
erqdaiiied  by  an  irKreaae  of  the  signal  to  noise  ratio  (SAfR)  for  small  onsets  due  to 
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beam  fonning  for  the  GERESS  dau  and  by  the  general  better  noise  conditions  at 
the  GERESS  site. 

At  first,  the  seasonal  pattern  was  investigated.  Fig.  3-2  shows,  how  the  de¬ 
tection  probability  of  GER^S  and  of  YKA  changed  with  respect  to  the  day  of 
year  frtra  OctobCT  18,  1991  to  February  11,  1993.  The  detection  capability  was 
averaged  for  10  days.  The  crosses  show  the  detection  capability  of  GERESS  and 
the  triangles  the  detection  capability  of  YKA.  The  two  arrays  differ  YKA  shows 
a  decrease  of  its  detection  c^idiility  of  about  25  %  for  the  summer  time,  viiien  the 
Great  Slave  Lake  is  open.  In  the  summer  and  early  autumn,  YKA  is  clearly 
less  sensitive  than  GERESS,  for  which  such  a  seastmal  effect  is  not  observable. 
The  decrease  of  the  d^ection  probability  of  GERESS  between  the  70'th  and  90*th 
day  of  the  year  was  caused  by  a  larger  outage  of  GERESS.  The  single  station 
KHC  situat^  in  the  same  geological  unit  as  GERESS  shows  a  similar  p^utem 
(Fig.  3-3). 

Bg.  3-4  shows  the  percentage  of  events  detected  by  the  three  seismic  obser¬ 
vatories  with  respect  to  epicaitral  distance.  Once  again,  triangles  are  for  YKA, 
crosses  for  GER^S,  and  circles  for  KHC.  The  decrease  for  GERESS  and  KHC 
for  events  between  60°  and  70°  epicentral  distance  around  the  North  Pacific  sub- 
duction  zones  is  due  to  die  large  number  of  smaller  events  (mi,  <  4J).  These 
events  were  detected  by  local  and  regional  stations  and  NEIC  could  determine  a 
/»£,  value.  Clearly  seen  is  the  decrease  m  and  around  die  Earth's  shadow,  but  also 
the  relatively  high  detection  capability  for  PKP  phases  between  I3S°  and  I5S° 
epicentral  distance.  Surpriskig.  the  highest  detection  capalnHty  in  this  plot  can  be 
found  for  PKP  (xisets  at  the  GERESS  site.  This  partly  reflects  the  low  network 
detection  ca^iability  for  smaller  events  (m^  <  4S)  in  the  South  West  Pacific.  It 
also  demonsrates,  that  sensitive  seismic  stations  in  epicentral  distances  around  the 
PKP  caustic  can  improve  the  detection  capability  drastically.  GERESS  rqxnted 
many  PKP  phases  in  the  investigated  time  range ,  which  cotild  not  be  associated 
to  any  event  by  NEIC 

At  next,  the  detection  capabUiOr  of  CARESS  during  the  time  of  day  was 
investigated.  Fig.  3-5  shows  die  detectimi  probability  of  GERESS  (upper  line) 
and  of  KHC  (lower  line)  widi  reflect  to  lo^  time  for  the  investigated  time  pe¬ 
riod.  A  decrease  of  the  detection  capability  at  daylight  hours  can  be  observed  at 
bodi  stations.  This  is  the  influence  ctf  an  increase  of  man  made  noise  and/or  of  the 
local  and  r^onal  artificial  seismicity  which  masks  smaller  teleseismic  onsets. 
This  interpretation  of  the  obsoved  eifea  is  confimied  by  Bg.  3-6,  vdiere  the 
detection  probability  during  the  whote  day  is  shown  only  for  wt^dng  days.  Now 
a  decrease  of  the  detection  probidrility  of  about  6  %  for  GERESS  and  about  5  % 
for  KHC  at  dayli^  hours  can  be  obs^ed. 
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3^  GERESS  mislocation  vectors 


GERESS  as  a  small  aperture  array  (~  4  km)  has  only  a  limited  resolution  for 
estimating  slowness  and  backazimuth  of  teleseismic  onsets.  After  association  of 
the  reported  onsets  to  located  events  by  NEIC,  the  reported  values  which  include 
slowness  ip)  and  backazimuth  {BAZ),  can  be  compared  with  the  theoretical  p  and 
BAZ  values.  NEIC  locations  were  used  to  calculate  the  epicmtral  distance  and  the 
theoretical  BAZ.  The  LASP91  travel  time  tables  (Kennett  and  Engdahl,  1991)  were 
used  to  calculate  the  theoretical  slowness  of  the  first  arrival  (P,  PKP)  from 
the  theoretical  qncentral  distance  and  the  reported  depth  (NEIC).  To  reduce  the 
scatter  all  p  and  BAZ  values  of  first  arrivals  with  a  travel  time  residual  /  tpbs  * 
hheo  /  ^  ^  ^  ^  slowness  mislocation  vector  in  the  slowness  space  of 

more  than  ±  4.0  secl°  were  omitted.  With  the  first  arrivals  investigated  during  the 
GSETT<2  period  (Sdiweitzer,  1992)  and  the  onsets  reported  between  Octobn- 18, 
1991  and  February  18,  1993,  a  total  of  3330  slowness  and  backazimuth  values 
could  be  compared  with  their  thec^etic^  values.  The  theoetical  coverage  of  the 
slowness  ^aoe  of  these  arrivals  can  be  seen  in  Fig.  3-7.  This  figure  should  be 
oonquaed  widi  3-1  to  locate  well  observed  source  regioQS  in  the  sbwness 
space. 

Hg.  3-8  shows  the  observed  mislocation  vectors  in  the  slowness  space.  The 
ini^k)cation  vectors  have  been  averaged  for  bins  with  a  width  for  /?  of  /  seel"*  and 
a  width  for  BAZ  of  10^.  In  the  central  part,  the  width  of  the  bins  for  BAZ  was  in¬ 
creased  to  20°  and  30°.  The  theoretical  value  is  always  indicated  by  a  circle.  The 
mean  slowness  residual  of  all  observations  is  s  -0.45  secl°  with  a  standard 
deviation  Ofjp  «  1.15  secl°  and  the  mean  BAZ  residual  is  dBAZ  -  4.76°  with  a 
standard  deviation  of  %  22.80°.  Because  of  the  large  standard  deviatimis, 

the  data  have  to  be  reduced  to  the  most  significant  observations. 

Bg  3-9  now  shows  only  mislocation  vectors,  which  were  estimated  from  at 
least  5  observations,  ft  is  obvious  that  the  remaming  mislocatioii  vecton  show 
systematic  effects.  The  BAZ  residuals  change  bt^een  positive  residuals  observed 
from  ~  350°  to  -  30“  and  from  -  90“  to  /  70“  and  negative  residuals  from  -  40° 
to  -  80°  and  from  -  180°  to  -  340°.  The  observed  slowness  valttos  are 
^stematically  too  small  widi  reqject  to  1ASP91.  At  the  moment  is  not  clesr,  how 
far  the  small  p  values  are  an  effect  of  the  Earth’s  structure  or  an  attract  of  dw 
implememed  fk-analysts,  the  configuration,  the  local  geaiogy,  and  die  to- 
pogrsidiy  at  the  array  site. 

Since  Jaiuary  1992,  YKA  and  GERESS  exchange  rqgularly  their  phase 
readings  via  dectronic  nudL  Therdbie,  the  reported  p  and  BAZ  vahies  YKA 
are  available  at  Bochum  and  we  can  compare  them  widi  their  theoretical  values 
respectivdy.  Fig.  3-10  shows  the  YKA  ndtlocstion  vectors  of  3679  onsets  for  die 
dme  period  between  Jamuuy  6, 1992  and  February  18, 19^.  The  mean  slowness 
error  YKA  is  Oj02  jec/“  with  a  standard  devittion  oi  Of^  »  0.46  seer*  and  die 
man  BAZ  error  is  0.6/“  widi  a  standard  devistion  of  OfQAZ-  53/“- Cotaraiy  to 
G9SRESS,  YKA  had  been  calibrated  and  therefore  only  oorreexed  p  and  BAZ 
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values  were  reported  and  its  remaining  mislocation  vectOTS  are  much  smaller  than 
for  CERESS.  fo  addition,  the  resolution  of  the  larger  YKA  array  is  much  better. 

Most  of  the  mislocation  vectors  of  GERESS  in  Rg.  3-9  are  small  with  re- 
speci  to  their  standard  deviation.  Therefore  it  is  not  useful,  to  correct  the  t^iserved 
p  and  BAZ  values  with  these  mislocation  vecton.  Fig.  3-11  only  shows  the 
mislocation  vectcns  for  which  die  theoretical  value  is  outside  the  error  ellipses 
calculated  with  die  standard  deviations  of  these  vectors  respectively.  At  the  mo¬ 
ment,  the  number  of  the  remaining  mislocation  vectors  with  "significant"  devia¬ 
tions  from  the  theoretical  p  and  BAZ  values  is  too  small,  to  use  them  routinely  in 
the  daily  reanaiysis  work.  Until  a  more  complete  set  of  mislocation  vectors  be¬ 
comes  available,  the  mean  standad  deviations  (o^  -  1.15  seer  and  > 

22B0°)  can  be  used  to  wei^  the  GERESS  slowness  and  bkdeazimuth  values  in 
location  routines,  which  use  this  additional  infoimation  (e.  g.  IMS). 


3.4  Travel  time  and  amplitude  residuals  of  GERESS 

In  this  diapter,  die  observed  bias  of  absolute  onset  times  and  amplitudes  at 
GERESS  are  diaaiascxl.  To  calculate  a  mean  travel  time  residual  for  GERESS,  All 
events  for  adiicfa  a  GERESS  arrival  significantly  contributed  to  the  location  must 
be  omitted.  Therefore  only  onsets  of  such  events  (A  2  20°)  were  used,  a^ikh  had 
fftf,  ^  5.0  (NEIC)  and  wiudi  were  located  widi  more  than  J9  dining  P  phases. 
Additionally,  all  onsets  with  an  absolute  value  of  the  travel  time  residual  larger 
than  S.0  sec  were  rejected.  Then  the  mean  travel  time  rendual  of  1235  GERESS 
(diases  was  estimated  at  ^  ^  >  standard  deviation  of  ~ 

sec.  Several  effects  can  eqilain  this  bias: 

The  sensitivity  of  a  station  is  die  most  important  contribution  to  the  travel 
time  residual  which  »  demoostialed  widi  Fig.  3-12:  the  mean  bavd  time  residuals 
of  other  arrays  and  some  European  stations  around  GERESS  are  shown.  Residuals 
are  plotted  widi  reflect  to  the  percentage  of  observed  EEMK  events.  The 
comparison  is  based  on  the  oontril^on  of  seismk  stations  in  the  time  period 
from  April  ISKX)  to  Fdiraaiy  18,  1993.  The  calailsled  percentage  oi  obaerved 
events  gives  a  rough  idea  of  the  sensitivity  of  die  stations.  All  esthnated  parame¬ 
ters  are  lisled  in  Table  3-1.  ft  is  obvious,  that  all  tafpi  sensitive  anrqv  haw  nega¬ 
tive  travel  time  residuals  and  dut  the  single  stations  show  a  sim^  trend:  die 
more  eensitive  die  more  negative  is  the  statton's  travel  time  bias.  This  obsenr«ion 
is  confirmed  by  Grand  (19^),  who  showed  that  the  travel  time  residuals  of  dif¬ 
ferent  stations  in  die  KC  bulletins  ocmeltfe  with  die  anpHfication  of  die  statkm. 
Thb  rerah  leads  to  die  concittnoti,  that  eventually  die  sooioe  time  of  larger  events 
is  systenmically  biased  to  later  dmes  due  to  die  large  number  d  less  sensitive 
Uatkxis  which  observed  diese  events. 

Phtally  the  »wg»«**“**  bias  of  (^RESS  will  be  disenised  ,  The  same  events 
as  in  die  tnv^  tinie  residual  study  have  been  used  to  calculate  die  mean  mnpii- 
tude  re^dual  for  CSRESS.  Due  to  the  huge  nnndier  of  PKP  onaets,  the  nundier 
mt,  vatuei  is  only  794.  For  dieie  angifitudes,  a  mean  res  of 
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dard  deviation  of  =  034  was  estimated.  This  relatively  large  magnitude 

bias  camoc  easily  be  explained.  At  the  moment,  we  assume  that  the  GERESS  mfy 
residual  is  a  summation  of  different  effects:  GERESS  has  its  best  SNR  for  P  on¬ 
sets  in  die  frequency  raige  between  IjO  Hz  and  3.0  Hz  (lost,  1992).  Therefore, 
this  narrow  frequency  band  is  often  used  to  analyse  a  teleseismic  arrival.  It  is  un- 
portam  to  note  that  noany  of  diese  onsets  would  not  be  vmbte  on  a  WWSSN  short 
period  instrument  (mean  frequency  ~  1.0  Hz)  at  the  GERESS  she.  This  narrow 
band  fUtering  increases  the  SI^  but  decreases  the  mean  amplitiidrs,  Additionatty, 
some  signal  loss  can  occur  due  to  die  beam  forming.  Ehally  at  die  GERESS  key 
station  GEC2,  we  do  not  observe  any  amplification  dire  to  a  sedimentary  l^er. 
All  these  points  tend  to  decrease  the  amplimdes  observed  at  the  GERESS  site.  We 
assume  as  die  most  impoitant  effect  that  GERESS  ami^itudes  are  measured  at 
relative  short  periods  wi^  reflect  to  most  other  stations. 

Using  log  A/T  instead  of  log  A  in  the  defirndoo  of  should  compensate 
fOT  different  instrumentation  and  observmg  periods.  But  the  mi,  residuals  of  all 
statians  used  in  die  travel  time  study  suggest  that  this  compensation  is  not  suf- 
fickm.  Fig.  3-13  shows  die  reskhials  of  different  stations  plotted  with  respect 
to  the  mean  reported  period  T  to  crxnpute  )og  A/T  (all  parameien  are  given  in 
TdUe  3-1).  The  stations  with  periods  Im  than  1.0  sec  ihaw  a  n^adve  reaid- 
uaL  Signal  loss  due  to  beam  forming  and  reading  parameters  from  a  beam  trace 
can  have  some  addhiotuil  effect  so  that  most  smy*  have  a  negative  lesiduaL 
All  diese  dtfa  are  influenced  by  the  local  geology,  but  if  only  stations  are  com¬ 
pared  wtudi  are  in  the  same  geological  unit,  the  main  influeinoe  seenis  to  be  the 
period  at  vdiidi  amplinides  are  measured  (e.  g.  GEC2,  KHC,  WET,  NMS,  FRU, 
MOX,VRAC). 

Condiisioii 

In  this  rqxirt  teteseismic  phases  detected  at  die  GERESS  array  were  evalu¬ 
ated  widi  leqrect  to  the  EI>R  bulletins.  No  seasonal  effect  of  dre  detection  capa- 
bflity  for  GERESS  was  found.  As  expected  dre  detection  capability  of  (XRESS 
duBiges  widi  epioeittral  ditianor,  for  events  in  the  Earth's  core  shadow  between 
100^  and  12(y*  qpioentral  distance  CARESS  observed  only  10  %  to  20  %  of  sH 
coirespondiiig  QMR  events.  In  contrast,  for  PKP  arrivals  GERESS  observed 
**  an  EIXt  events  Are  to  the  PKP  caustic  effect  The  daily  variatkn 

('«  d  %)  of  the  detection  capability  for  GERESS  could  be  associared  with  an  in- 
crease  of  cuitutal  and  iiidua&ial  nmse  during  working  hours. 

AldKNi^  CnSRESS  widi  its  smaU  aperture  has  only  a  limited  reatdntion  in 
mraairint  slownen  and  barkinimith  of  tekarejanic  onsets,  ^stematic  miikicar 
tkm  vectors  were  found  fOT  GERESS. 

Finally,  travd  time  and  amplitude  residuals  of  C^BRESS  were  discuned. 
The  iiegadvetravd  dine  residual  of  CERESS  can  inorty  be  exptmed  by  dre  hatfi 
aearitivitytrfdre  atn^.  The  n^;ative  residual  of  GERESS  aeema  to  be  a  aum- 
mataon  of  different  effecre  which  aE  to  dimmiah  ampHnidei  with  leapeet  to 
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other  stations.  For  both  -  travel  tizms  and  amplitudes  •  the  GERESS  residuals  are 
conyarable  with  residuals  hom  other  hi^  sensitive  arrays. 
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Table  3<1:  Estimated  travel  time  residuals  (t,^  and  magnitude  residuals  (mb  fes) 
for  diffensnt  arrays  and  European  stations  with  the  luin^jer  of  observations  (#  of 
events,  #  of  015,  T).  Also  list^  are  the  percentage  of  observed  events  whh  respect 
to  the  total  number  of  events  in  the  investigated  time  period  (%  of  events)  «id  the 
mean  poiod  of  reported  an^litude  measurements  to  calculate  0%  (mean  T). 
Additionally  given  are  the  standard  deviations  (Oj;  Of^bi  crp)  respectively. 


Station  t|^  %of  #of  mb^s  or|nb  meanT  crp 

events  events  [sec]  nib,T 


essB 


0.73  110.1 


11.4 


BRC 

1.01 

022 

472 

1933 

-0.03 

022 

1.36 

0.37 

CLL 

-0.11 

023 

51.7 

2104 

0.04 

024 

1.45 

0.38 

EKA 

-0.01 

022 

31.9 

r0.14 

027 

1.06 

0.39 

ENN 

027 

0.33 

28.0 

1139 

0.04 

027 

1.01 

0.24 

FUR 


GBA 


GEC2 


GRF 


KAF 
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KHC 


KSP 


UU 


LOR 


0.54 


0.40 


-OJO 


0.45 


-1.06 


•0.84 


0.08 


0.24 


0.61  19.6 


0.43  16.6 


0.44  33.9 


0.80  70.7 


0.47  42.3 


1.06  45.0 


0.85  45.1 


0.24  24.3 


0.30  51.1 


0:22  32.7 


TNS 


VKA 
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WET 
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0.16 

0.02 

0.19 
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Bgore  3-5:  The  probability  to  observe  an  event  (A  >  2(P)  during  the  time  of  day 
in  the  investigated  time  period  (October  18,  1991  -  Fdnuary  11,  1993)  for 
C£R1^  (upper  line)  and  KHC  (lower  line). 


Rgna  3-6:  Aa  Fig.  3-5,  bat  only  for  woricn^  days  in  Germany. 
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FIgiife  3>8:  OERESS  mean  nualocation  vectois.  The  symbols  are  plotted  at  die 
ttMoradeal  vdue.  Fbr  ftaidier  detail 


Hfure  3-10:  As  Bg.  3-9,  but  for  YKA.  Giveo  is  die  numbers  of  observations  for 
the  vector  and  the  ttndaid  deviatioos  are  shown  as  dlipae. 


MEAN  STATION  RESIDUAL  [SEC] 

Figure  3-12:  Percentage  of  NEIC  events  observed  at  dilTerem  seismic  staticms 
venus  the  station  travel  time  residuals  req)ectively.  For  more  details  see  text 
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Figure  3<13:  Mean  reriduala  of  different  sttfions  fribtied  widi  re^Msct  to  die 
mean  measured  periods. 
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4.  STEEP  ANGLE  REFLECTIONS  FROM  THE 
EARTH’S  CORE  OBSERVED  WITH  THE  GERESS 

ARRAY 

Johannes  Schweitzer 


4.1  Introduction 

The  systematic  search  for  unusual  phases  has  found  more  interest  in  the  last 
years.  Hiese  (4»ervations  can  be  used  either  to  study  the  spherical  structure  of  the 
Earth  or  serve  as  an  additional  informatitm  in  tomogra;^c  studies.  Detailed  in¬ 
vestigations  of  rqxffted  f^ases  to  the  ISC  and  NEIC  (Doombos  et  aL,  1991; 
Hwang  and  Qayum.  1991)  diow,  that  many  later  arrivais  in  seismograim  can  be 
identified  and  a^ysed. 

The  scientific  goal  of  the  International  Seismological  Observing  Period 
(ISCX*)  is  to  collect  not  only  observations  of  hi^  quality  of  the  standard  seismic 
phases  but  also  to  complete  those  parts  of  travel  time  curves,  which  are  rx>t  yet 
confirmed  with  sufficient  observati<»is.  Because  these  [rfuBes  are  often  small  and 
disturbed  by  odier  more  domirumt  onsets,  i^iase  identification  and  interpretation 
has  to  be  d^  very  carefully.  One  amy  to  distinguish  between  different  onsets  is 
the  usage  of  array  tedmiques. 

Fig.  4-1  shows  the  theoretical  travel  time  curves  for  a  surface  focus 
(IASP91  model,  Kennea  and  Engdahl,  1991)  and  in  Fig.  4-2  travel  times  repotted 
to  the  ISC  and  accumulated  for  10  years  are  shown  (Doombos  et  aL,  1991).  The 
core  reflections  PcP  and  SeP  are  usually  tx>t  observed  or  not  repotted  at  epwentral 
distances  smaller  than  20°  and  for  PKiKP  at  epicentrai  distances  smaller  than 
lOff*.  Id  dus  study  the  most  sensitive  seismic  station  in  Ceiural  Europe,  die 
GERESS  array,  is  used  to  search  for  such  steep  angle  reflections  from  the  Earth's 
core. 


4J2  PcP  observatioiis  in  less  tfian  3S^  epicentrai  distance 

In  a  seardi  for  steep  angle  PcP  phases,  the  GERESS  data  base  from  Afnil 
1990  to  Fdmiary  1991  was  analysed.  The  selectit^  criteria  were  a  value  o( 
larger  than  5.0  (NEIC)  and  an  epicentrai  distance  less  than  35°.  For  this  dK>sen 
time  period,  25  events  were  found  in  the  GERESS  data  base.  The  source  parame¬ 
ters  ate  summarised  in  Table  4-1  and  the  locations  of  the  events  ate  shown  in  Fig. 
4-3.  Only  in  fhree  cases  (filled  triangles),  a  clear  PcP  onset  was  found,  bt  one  esse 
a  small  onset  could  be  observed  but  this  onset  could  not  be  identified 
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T^le  4-1:  Eveitts  used  in  this  PcP  study.  All  parameters  for  events  in  1990  are 
from  the  ISC  bulletins,  for  the  1991  event  from  NEIC. 


DATE 

TIME 

LAT 

n 

LON 

n 

DEPT 

H[km] 

mb 

A 

PcP 

1990  Apr  20 

23:30K)5.1 

40.12  N 

22 

5.0 

20.7 

im 

1990  May  05 

07:21:29.6 

40.74  N 

15.74  E 

9 

5.1 

83 

— 

B 

1990  May  30 

10:40^)6.4 

45.85  N 

26.66  E 

89 

6.4 

9.3 

B 

1990  May  31 

11:17:48.2 

45.81  N 

26.77  E 

90 

5.9 

la 

1990  Jun  09 

17:24:52.6 

40.16  N 

29.71  W 

10 

CB 

313 

B 

1990  Jun  16 

02:16:21.1 

39J7N 

2035  E 

32 

53 

10.8 

m 

09K)2:14.8 

36.65  N 

49.83  E 

11 

5.1 

28.9 

yes 

1990  Jun  21 

12:17:32.2 

36.94  N 

49.48  E 

41 

5.8 

283 

B 

1990  Jul  06 

19:34:54.3 

36.91  N 

49.30  E 

45 

53 

28.4 

1990  Jul  09 

11-.22:18.1 

34.94  N 

26.63  E 

41 

5.0 

16.9 

B 

1990  Jul  24 

17:59-.50.7 

32.84  N 

4832  E 

24 

5.3 

303 

B 

1990  Sep  09 

5639  N 

34.39  W 

19 

5.1 

29.6 

B 

1990  Sep  09 

06:35:22.7 

56.48  N 

34.44  W 

10 

m 

29.6 

B 

12:17:37.0 

70.75  N 

13.80  W 

10 

5.0 

253 

B 

1990  Sep  15 

23K)7:43.2 

64.75  N 

17.50  W 

18 

5.3 

23.0 

B 

1990Oct  08 

23:42:42.5 

1932  E 

10 

a 

6.1 

B 

1990Oct  24 

14:57:58.2 

7336  N 

54.67  E 

0 

m 

30.4 

yes 

1990  Nov  12 

16K)0:49.6 

40.31  N 

51.63  E 

44 

53 

28.1 

— 

1990  Nov  27 

04:37:57.0 

43.87  N 

16.67  E 

19 

53 

5,4 

^31 

1990  Nov  27 

04:51:36.5 

43.90  N 

16.65  E 

10 

5.0 

5.4 

B 

1990  Dec  13 

00-,24a5.7 

37.31  N 

15.43  E 

53 

11.6 

1990  Dec  14 

03:21:26.9 

39.35  N 

275 

4.9 

9.6 

yes 

1990  Dec  16 

15:45:40.3 

4137  N 

43.72  E 

28 

5.3 

22.4 

B 

1990  Dec  29 

15:34:17.2 

3239  N 

47.86  E 

34 

5.1 

303 

1991  Feb  26 

07-.25:47.2 

40.17  N 
_ 

13.82  E 

401 

U 

8.7 

B 

obviously  as  PcP.  For  all  other  evoits  no  PcP  ons^  were  found.  This  illustrates 
again  the  difGculties  to  observe  PcP  under  steep  angle  conditions. 

The  most  prominent  PcP  observaticm  was  from  the  nuclear  test  at  Novaya 
Zemlya  cm  October  24,  1990.  The  seismograms  (Fig.  4<4)  are  tkxninaied  by  ap¬ 
proximately  three  minutes  of  large  amplitude  after  the  first  onset  and  an  addi¬ 
tional  sharp  (xis^  after  about  200  seconds.  The  observed  P  coda  is  mostly  due  to 
the  mantle  triplication's,  P  energy  reflected  multiply  at  the  Earth's  surface  (e.  g. 
PP,  pPP,  sFP,  VPPf ...),  and  energy  scattered  at  lateral  h^erogeneities  in  u|^)er 
mantle  ^  crust  The  sharp  onset  is  the  PcP  phase,  whidt  can  be  identified  even 
on  single  unfilteicd  traces  of  the  GERESS  array.  Fig.  4-5  exhibits  PdP  in  more 
detail  with  the  GERESS  PcP  beam  trace  on  the  top.  PtP  has  the  classical  simple 
pulse  fonn  as  expected  from  an  erqrlosion  source  contrary  to  the  complex  seis¬ 
mogram  of  the  P  onset  (Fig.  4-6).  Because  of  the  simple  PcP  pulse  fonn  the 
complex  P  pulse  form  can  only  be  explained  with  hderogeneities  akxtg  die  P  ray 
path  from  Novaya  Zemlya  to  GERESS. . 

Fig.  4-7  shows  die  seismograms  of  an  event  in  Western  Iran  (June  21, 1990) 
with  Mtnilar  strength  and  qiioentral  distance  as  the  last  nuclear  test  of  the  Soviet 
Union.  The  seismograms  again  diow  a  clear  PcP  onset  on  all  GERESS  single 
traces  but  the  P  coda  is  much  smaller.  However,  additionally  clear  depth  j^hase 
energy  (i^,  sP)  can  be  seen  (source  depth  41  km).  The  energy  increase  after  PcP 
can  be  explained  as  pPdP  and  sPdP  energy.  Fig.  4-8  shows  die  Butterwoidi  band 
pass  (0.5  Hz  -  4.0  Hz)  Altered  GERESS  beams  for  P  and  PcP. 

Some  days  later  (July  6,  1990)  another  event  occurred  in  the  same  source 
region.  The  seismograms  (F^  4-9)  show  -  as  in  the  case  of  the  nuclear  test  -  a 
larger  P  phase  coda  and  significant  eneigy  for  the  surface  reflections.  Once  again, 
an  onset  at  the  presumed  PcP  arrival  time  can  be  observed  at  the  GERESS  array. 
But  the  fk-analysis  of  this  ons^  yields  a  phase  velocity  whidi  is  to  low  to  be 
<»cplaiiied  as  a  ptuae  reflected  from  the  Earth's  core.  Tids  exangile  demonstrates 
again  the  advantage  to  use  an  array  to  identify  arrivals.  The  observed  difference 
between  the  Pd*  observations  fitxn  two  events  in  the  same  source  r^ion  can  only 
be  explained  by  different  focal  mechanisms  of  the  two  events.  Fig.  4-lO.a  shows 
the  Harvard  CMT  solution  (from  the  ISC  bulletin)  of  die  first  event  and  Fig.  4- 
10.b  shows  the  Harvard  solution  of  the  second  event,  reflectively.  The 
second  event  has  a  very  small  Fd*  radiation  and  the  large  P  phase  raditt^  (T- 
axes)  in  direction  to  GERESS.  Bodi  facts  reduced  the  possibility  to  obsove  PcP. 

It  is  not  surpiising,  that  Pd*  is  difficult  to  observe  in  the  cross  over  distance 
widi  die  S  fdiase  in  the  vicinity  of  20°  epicentrai  distance.  The  PdP  reflectkm  co- 
efAcient  (Fig.  4-11)  at  the  core  mantk  boundary  (C3llB)  deceases  whh  deoeas- 
ing  epicentrai  distance  and  the  S  energy  of  larger  evoits  noimally  dominates  the 
seismogram  on  an  3  components  at  these  diatances.  Another  possibility  to  observe 
Pd*  «tists  at  distances  up  to  about  1(P  epicentrai  distance,  this  case,  the  array 
tecfankpie  can  be  used  to  s«9ipress  the  S  coda  and  surface  waves  by  vdodty  filter¬ 
ing.  The  excitation  of  surface  waves  is  normally  small  for  intermediate  <»r  deep 
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events.  In  this  study  four  of  such  events  at  tegional  distances  were  investigated 
and  tite  data  of  two  events  are  discussed  in  the  following. 

The  first  event  was  the  large  Romanian  earthquake  occurring  on  May  31, 
1990  (A  ~  9.4^*,  Hq  >  90).  The  focal  medianism  of  this  event  (Hg.  4'12ui)  shows 
a  good  PcP  radiation  but  also  good  P  and  S  phase  radiation  in  the  direction  of  the 
GERESS  array.  Gmsequently  the  observ^  seismograms  show  ixKxe  or  less 
constant  amplitudes  over  several  minutes  (Fig.  4<13).  Therefore  it  is  very  difficult 
to  identify  the  PcP  onset 

The  last  example  in  this  PcP  study  is  the  event  of  December  14,  1990  (A  = 
9.6’^,  Hq  ~  275  km).  The  focal  mechanism  (Rg.  4-12.b)  cornpiled  from  dK  first 
mtHion  data  in  the  ISC  bulletin  shows  clearly,  that  the  P  and  S  radiatitm  in  the 
direction  to  GERESS  {azimuth  =;  353°)  was  relatively  smalL  Hus  is  confinned  by 
the  GERESS  seismograms  (Fig.  4-14).  After  the  P  ccxia,  no  significant  S  phase  is 
seen  and  Lg  waves  oi  otho*  surface  waves  are  missing,  hi  the  theoretical  time 
window  of  PcP  a  small  arrival  can  be  seen  in  the  unfiltered  seismograms  (Fig.  4> 
IS),  which  was  also  detected  by  the  automatic  GERESS  data  analysis,  l^th  a  fk> 
analysis  a  ray  parameter  p  ~  1.9  sed°  and  a  backazimuth  {BAZ)  of  118.8°  could 
be  estimafed  ^^thin  the  resolution  of  die  GERESS  array  for  hi^ier  ve- 
locities,  these  values  are  in  agreement  with  the  theoretical  values  (p  s  0.9  sec/°, 
BAZ  s  1723°).  Approximately  5  sec  later  a  PP  onset  of  an  evem  in  Taiwan  (mif 
~  52)  can  be  expected  at  die  GERESS  array  (theoretical  values  of  p  ss  82  sec/° 
and  BAZ  s  613°).  With  the  result  of  the  fk>analysis  this  hypothesis  for  the  ob- 
served  onset  can  be  rejected  and  die  ons^  can  be  cleariy  identified  as  a  PcP  ob- 
sovation  in  ,4  «  9.6°. 


43  ScP  ObservBtHMis  in  less  than  20°  epicentral  distance 

If  the  SV  radiation  of  events  in  direction  of  the  core  is  large  enough,  the  S  V 
to  P  converted  reflection  fhnn  the  CMB  (ScP)  can  be  trfiserved.  The  arrival  time 
of  ScP  is  approximately  200  sec  after  PcP  and  therefore  less  disturbed  by  S 
idiases  or  surface  waves.  Additionally  the  reflection  coefficient  for  ScP  at  die 
CMB  is  larger  as  in  the  case  of  PcP  (Fig.  4>I6),  in  die  observed  cfuceotral  dis¬ 
tance  rmge  from  12°  to  18°  the  reflection  coefficient  of  ScP  is  coixqiardile  with 
die  PcP  reflection  coefficient  between  23°  and  35°  epicentral  distance.  W^iout 
any  systematic  search  for  ScP  in  the  GERESS  data  base  some  ScP  observations 
can  be  shown  in  dus  lepoiL  The  GERESS  detector  found  onsets  at  die  theoretical 
arrival  time  of  ScP  for  three  recent  events  in  qncentral  distances  less  than  20° 
(see  TaUe  42).  The  results  of  the  fk-analysis  show,  that  these  arrivals  have  the  ■ 
light  BAZ  and  slowness  as  expected  for  ScP.  Rg.  4-17  shows  the  Buttowmth 
band  pass  fiheied  traces  {0.7  Hz  -  3J)  Hz)  oi  diese  observations,  once  for  die 
C^ERESS  beam  and  once  for  the  data  of  the  GERESS  reference  site  GEC2. 
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Table  42:  The  events  used  for  the  ScP  study.  The  parameters  are  from  NEIC 
weekly  listings. 


DATE 

TIME 

LAT  n 

■BH 

MBM 

mh 

li 

1992  Nov  18 

21:10:40.9 

38.30  N 

22.43  E 

10 

5.9 

12.3 

1992  Nov  21 

05:07:22.8 

35.93  N 

22.44  E 

70 

5.9 

14.4 

1993  Jan  09 

24:52:44.5 

35.88  N 

29.94  E 

32 

17.6 

4.4  A  PKiKP  observation  with  the  GERESS  array  in  52.4® 

PKiKP  is  rarely  observed  at  epicentral  distances  less  than  90°.  Due  to  the 
«nall  reflection  coefficient  (Fig.  4-18)  of  PKiKP  at  the  iimer  core  bcHindary  (ICB) 
of  less  than  5  %,  only  small  amplitudes  can  be  expected  for  such  steep  angle  re- 
flecticxis.  especially  since  the  wave  also  has  to  penetrate  additionally  twice  the 
CMB.  The  large  Qiinese  nuclear  test  (mi,  =  65)  <ni  May  21,  1992  at  05K)0  UT 
allowed  the  analysis  of  this  phase.  Fig.  4-19  shows  that  the  S  phase  radiation  of 
this  etqplosion  was  very  small  (as  erqsected)  and  that  the  P  radiation  to  the  Earth's 
core  was  relatively  good  -  as  the  clear  PcP  ons^  shows.  The  search  for  PKiKP 
onsets  for  this  event  at  the  Scandinavian  arrays  (ARCESS,  FINESA,  NORESS), 
at  the  Grafenberg  array  (GRF),  and  at  seven  broad  band  stations  in  Central 
Eurt^  (BFO,  BUG,  CLZ,  FUR,  KSP,  TNS,  WET)  was  not  successful,  only  at 
the  GERESS  array  a  PKiKP  onset  could  be  identified 

The  GERESS  short  period  vertical  traces  were  filtered  with  a  Butterworth 
band  pass  niter  (1.0  Hz  -  3.0  Hz)  since  the  GERESS  array  has  the  best  signal  to 
noise  ratio  (SNR)  for  P  onsets  (lost,  1992)  in  this  frequency  band.  On  the  single 
traces  of  the  GERESS  array,  no  clear  onset  is  vitible.  But  tte  fk-analysis  (Fig.  4- 
20)  of  the  PKiKP  expected  time  window  shows  that  at  this  time,  energy  arrived 
wiA  the  high  phase  velocity  of  105.9  Msec  (1.05  secl°)  from  the  East  (BAZ  = 
80.9°).  The  diffrxence  to  the  theoretical  values  (103.0  Msec  or  p  =  1.09  secf°  for 
the  IASP91  travel  time  tables,  BAZ  ^  68.0°)  b  surprisingly  small  with  respect  to 
the  accuracy  of  slowness  measurements  of  the  GIQIESS  array  (see  Chapter  5  in 
tfab  repmt).  Additknal  investigations  widi  vespagrams  and  a  moving  window-fk 
analysb  (Gestermatm,  1992)  to  estimate  slowness  and  azimuth  of  thb  oroet  con¬ 
firm^  the  former  result:  very  high  velocities  and  a  BAZ  from  the  EasL  Hg.  4-21 
shows  the  GERESS  beam  calculated  with  the  results  of  the  fk-analysb  (Fig.  4-20) 
for  the  PKiKP  onset  in  the  middle.  At  the  top,  the  SNR-tzdce  (calculated  as  rela¬ 
tion  between  the  short  time  (1  sec)  to  the  long  time  (30  sec)  amplitude  average)  of 
the  beam  b  seen,  and  at  the  bottom  the  GERESS  reference  station  GEC2.  Addi¬ 
tionally,  tte  dieoretical  nCiKP  onset  time  of  the  1ASP91  travel  time  tables  is 
marked. 

A  PKiKP  observation  in  dib  distance  b  not  c«ily  interesting  because  of  the 
travel  time  tnfotnuaion  for  P  phases  travelling  thtm^  Ute  Earth's  core,  but  the 
existence  of  steep  angte  observations  of  PKilQ’  also  confirms  that  the  ICB  is  a 
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sharp  boundary.  The  amplitudes  of  PKiKP  can  be  used  to  investigate  the  reflec¬ 
tion  coefficient  at  the  ICB.  Shearer  and  Masters  (1990)  systematically  searched 
for  PKiKP  onsets  through  the  GDSN  data  base.  They  found  only  two  observati<»is 
of  this  rare  i^ase  m  a  time  period  of  S  years.  To  discuss  an^iitudes  of  PKiKP 
with  respect  to  the  reflection  coefficient  they  compiled  PKiKP  observations  from 
the  literature.  To  normalise  the  PKiKP  amplitudes,  it  is  best  to  give  the 
PKiKP/PcP  amplitude  ratio.  Therefore,  the  GERESS  PcP  beam  of  the  Chinese 
nuclear  test  was  processed  similaiiy  as  the  PKiKP  beam  (Fig.  4-22).  The  observed 
PKiKP  airplitude  at  GERESS  in  A  -  51.4^  is  %  of  the  PcP  amplitude.  This 
datum  has  beat  added  to  the  data  of  Shearer  and  Masters  (Fig.  4-23).  The  new 
observation  is  in  good  agreement  widi  their  data. 


4.5  Conciusitm 

The  seardi  for  steep  angle  reflections  from  the  Earth's  core  has  demon¬ 
strated  that  small  aperture  short  period  arrays  like  GERESS  can  not  only  be  used 
to  investigate  local  and  regional  seismicity.  Due  to  sites  with  hi^  SNR  ccmdi- 
tiotts,  diK  to  beam  forming  teduiiques  wiA  a  large  number  of  seismic  traces  in¬ 
cluding  some  redundant  informations,  and  due  to  fk-analysis  for  estimating  the 
phase  velocity  and  the  A4Z,  such  arrays  have  great  capabUittes  for  other  seis- 
mological  purposes.  The  capabilities  of  the  GERESS  array  can  be  also  used  to 
associate  later  arrivals  (e.  g.  PT’,  PKKP, ...)  to  the  apjnoptiate  events  by  unng  the 
BAZ  information  or  to  detect  small  evoits  at  telesdsmic  distances.  E^tecially  in 
connection  with  the  ISCNP  project,  the  advantages  of  small  aperture  arrays  like 
GERESS  can  be  used  to  sei^  for  rare  seismic  phases,  wfaidi  normaliy  cannot  be 
detected  or  cannot  be  identified  with  single  seismic  stations. 
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Figure  4-i:  Travel  time  curves  of  an  event  at  the  Earth’s  surface  Figure  4-2;  Reported  seismic  onsets  of  10  years  (1974  -  1984)  ISC 
(Kennett  and  Engdahl,  1991,  Fig  Al).  data  (Doombos  et  al.  1991,  Fig,  1)  for  shallow  events  (0  -  50  km). 


Flguxe  4-3:  Location  of  events  used  in  this  study  of  PcP.  For  events  with  filled 
triangles  GCRESS  observed  a  PcP  onset  The  event  with  the  filled  octaeder  shows 
a  questtonable  observation  of  PdP.  For  the  events  with  open  octaeders  no  PdP 
nfiyir  oonld  be  observed  at  QERESS.  The  map  is  an  aznnuthal  equidistant  projec¬ 
tion  with  CARESS  in  the  center.  Distances  with  respect  to  GERESS  are  given  for 
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Fig'uie  4-4:  Unnhercd  QERESS  record  of  P  phase  and  P  coda  group  from  the  last  Soviet  Union  nuclear  lest  at  Nnvaya  Zemlya  (October  24. 1900). 
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Figure  4-7:  Butterworth  band  pass  filtered  (OJz  -  2J  Hz)  seismograms  of  the  P, 
pP,  sP,  and  PdP  onsets  from  the  Western  Iran  event  (June  21, 1990). 
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Figure  4-9:  An  Fig.  4-7,  but  for  the  Western  Iran  event  on  July  6,  1990.  Note  the 
larger  P  ooda  with  respect  to  the  seismograim  in  Fig.  4-7. 
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E^ure  4-11:  Absolute  values  of  the  plane  wave  PcP  leflecticxi  coefficient  at  the 
Cl^  for  PREM  (Dziewonsid  and  Andenon,  1981)  as  fimction  of  the  angle  of 
inckkiioe  at  the  CMB.  Also  noted  are  the  epicentral  distances  ctntesponding  to 
the  ai^  of  incidence. 
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Figure  4-13:  Unfillered  seismograms  of  Ihe  Romanian  event  on  May  31, 1990  as  recorded  with  the  GliRESS  array. 


Figure  4-14:  As  Fig.  4-13,  but  for  the  Southern  Italy  event  «mi  December  14, 1990. 


Figure  4-16:  The  absolute  value  of  the  plane  wave  reflection  coefficient  of  SeP  (left 
y-oes)  IK  the  CMB  for  PREM  as  a  funaion  of  incidence.  The  broken  line  gives  the 
**•1*  etnergence  at  the  CMB  (right  y-axes)  and  the  corresponding  qncentral 
<&0HK9es  are  tnaaiced. 


Figure  4>18:  The  absolute  value  of  tin  plane  wave  reflection  coeflUuent  of  PfCiKP 
at  the  ICB  for  PREM.  The  incidence  angle  at  die  ICB  for  did  observed  PKiSP 
onset  at  CARESS  is  marked. 
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!992-142;05.17.00.000  +  4.00  SEC 
WIDE-BAND  SLOWNESS  ANALYSIS 
CONTOURS  IN  OB  BELOW  MAXIMUM  PEAK 
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Hgure  4>20:  The  resuh  of  the  flc^analysis  in  die  time  wincbw  of  the  PKiKP  onset 
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Rfuie  4-22:  At  Fig.  4-21,  but  for  the  PcP  onset. 
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Figure  4-23:  Comparison  of  the  PKiKP^*cP  observatiai  at  GERESS  with  other 
puHishfH  values  (modified  after  Shearer  and  Masters,  1990). 
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5.  GERESS  RECORDINGS  OF  PRESUMED  NUCLEAR 

TESTS  IN  1992 


Johannes  Schweitzer 


5.1  Introduction 

In  1992,  presumed  nuclear  explosions  were  detected  from  only  two  test  ^es: 
Nevada  and  Lop  Nor.  This  ch^er  is  the  contmuatioD  of  former  rqxMrts 
(Sdiweitser,  1991;  Sdtwehzer  et  d.,  1992)  on  nudear  tests  observed  widi  the 
GERESS  array.  The  known  or  presumed  nuclear  tests  in  1992  are  listed  in  the  fol¬ 
lowing  Table  5-1. 


Table  5-1 :  Presumed  or  confirmed  nuclear  tests  in  1992  and  the  number  of  obser¬ 
vations  at  the  GERESS  array  sorted  fix  the  diflFerent  test  sites. 


Country 

Test  Site 

E{»central 

Distance 

#  of  Tests 

China 

51.4« 

2 

2  1 

USA 

Nevada 

83.5® 

6 

mmnmm 

The  dififereooes  between  die  total  nundwr  of  tests  (8)  and  events  recorded 
with  a  dear  omet  (4)  is  caused  by  tests,  whidi  were  too  anall  to  produce  an  ob- 
servdde  s^^  at  die  C£RESS  she.  All  observed  GERESS  parameters  of  die  S 
midear  test  in  1992  are  listed  in  Table  5-2.  hi  the  following,  die  observed  data  wiU 
be  shown  test  site  by  test  site.  As  in  fijtmer  reports,  all  data  are  given  with  respect 
to  the  GERESS  key  station  C2EC2  and  all  parameters  are  measured  on  the 
(XRESS  beams. 


Table  5-2:  Observed  parameters  of  presumed  <m'  confirmed  nuclear  tests  as  recor¬ 
ded  at  GERESS  in  1992. 


Location 

Date 

Onset 

Time 

tsa 

Period 

fsec] 

SNR 

®b 

BAZn 

Nevada 

Mar  26 

16:42:29.0 

45.2 

1.29 

42.0 

5.6 

332.9 

3.53 

05:09:06.3 

334.8 

0.84 

364.0 

6.3 

65.8 

7.23 

Seo  18 

17:12:30.2 

1.4 

0.97 

5.3 

4.2 

348.3 

6.32 

15:16:30.8 

0.71 

1.8 

EE] 

(282.4) 

mam 

6.2 

1.33 

5.3 

4.4 

61.8 

7.691 

5.2  Lop  Nor 


In  1992,  China  continued  again  its  nuclear  test  program.  Two  presumed  nu¬ 
clear  tests  occurred  at  Lop  Nor,  which  both  were  observed  with  the  GERESS  ar¬ 
ray.  With  the  observation  of  May  26,  1990,  the  GERESS  data  base  now  contains 
three  exfdosioiis  from  the  Chin^  test  she.  The  first  explosion  in  1992  (on  May 
21)  was  the  largest  explosion  which  was  ever  observed  with  broad  band  stations 
wmid  wide.  This  initiated  frirfrier  sdmiological  investigations,  some  results  can  be 
found  in  Chapter  4  of  this  report. 

Fig.  S-1  shows  the  GERESS  BB13  3-component  registrations  of  the  May  21 
event.  The  horizontals  were  rotated  using  the  themetical  baclcazhnuth.  Clearly  seen 
is  the  laige  P  phase  amplitude  vduch  dominates  the  unfihmed  seismograms  (three 
top  traces).  Although  for  this  event  relatively  large  Love  waves  could  be  observed 
(see  the  three  WWSSN  s'nulated  traces  at  the  bottom),  the  classical  m3 /Ms  ratio 
(6.3  /  5.2  =  7.2)  at  the  GERESS  site  »ilL confirms  the  explorion  source.  To  don- 
onstrate  the  simple  pulse  form  of  reiativdy  hi^  frequency  for  this  large  event, 

5-2  shows  the  deconvolved  P  and  PcP  onsets.  In  b^  cases  a  long  periodic  down 
swing  after  the  first  onset  is  seen,  vrindi  can  be  addressed  as  an  e£^  of  the  tec¬ 
tonic  release  related  with  this  event  (Gao  and  Wallace,  1992).  Another  effect  of  the 
tectonic  rdease  seems  to  be  die  large  Love  wave  amplitudes. 

Fig.  5-3  shows  the  band  pass  filtered  data  of  the  second  Chinese  nuclear  test 
in  1992.  The  top  trace  is  again  the  best  GERESS  beam.  This  very  small  event  was 
detected  only  at  fisw  stations  in  Europe,  as  the  contributions  to  the  May  1992  EDR 
show. 

53  Nevada 

From  the  6  announced  nuclear  tests  at  the  Nevada  test  site  in  1992,  only  two 
were  detected  by  the  (XRESS  automatic  data  analysb  system.  The  4  other  tests 
were  too  small  to  get  a  detectable  onset  in  Central  Europe.  In  om  case  a  small  on¬ 
set  widi  a  a^nal  to  noise  ratio  (SNR)  Of  1.8  became  visiUe  <m  the  array  beam.  But 
it  is  very  dovAtfiil  to  claim  this  onset  as  an  observadoa  Therefore  the  observed  pa¬ 
rameters  of  this  onset  are  given  as  upper  limits  in  Tdile  5-2  and  the  parameteni 
calculating  the  beam  are  given  in  bradcecs. 

The  data  of  the  two  detected  Nevada  test  site  events  are  shown  in  Fig.  5-5 
and  Fig.  5-^.  The  uncertain  onset  of  the  small  explosi<m  of  Sqnendia'  23, 1992  can 
be  seen  in  Rg.  5-5.  Always  shown  are  the  band  pass  filtered  traces  of  the  CSRESS 
array  with  the  beam  trace  at  the  fop. 

After  3  years  of  momtoring  the  Nevada  test  site  with  the  GERESS  array,  the 
GERESS  data  base  now  contains  1 1  observations  in  the  magnitude  rai^  between 
W3  -  4.3$admif^  5.7. 
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Figure  5-4;  As  Fig.  5-3,  but  now  for  the  explosion  at  the  Nevada  test  site  on  March  26,  1992. 
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Figure  S-S:  As  Fig.  5-4.  but  now  for  the  nuclear  test  at  the  Nevada  test  site  on  September  18,  1992. 


6.  SYNTHESIS  OF  BROAD  BAND  DIS- 
PLACEMENT  SEISMOGRAMS  FROM 
BAND  LIMITED  DIGITAL  REGISTRA¬ 
TIONS 

Jan  W^ter 


6.1  Summary 

With  the  advent  of  high  resolution  analog/digital  converters,  it  is  no  longer 
necessary  to  use  broad  band  seismometers  in  order  to  achieve  complete  cov¬ 
erage  of  the  seismic  spectrum.  In  this  paper  a  method  is  presented,  with 
which  broad  band  records  of  true  ground  displacement  can  be  synthesized 
using  band  limited  digital  registrations  obtained  from  two  or  more  different 
standard  instruments  at  the  same  location.  This  process  involves  forward 
recursive  filtering  only  and  is  thus  suitable  for  on  line  real  time  application. 
Its  usefulness  is  demonstrated  in  an  application  to  the  instruments  at  the 
GERESS  reference  station  GEC2. 

6.2  Introducing  the  method 

At  the  beginning  of  the  1980s,  the  use  of  digital  broad  band  seismograph 
systems  was  strongly  advocated  [Harjes,  1981]  to  achieve  both  international 
standardization  of  seismological  records  and  equal  coverage  of  the  whole  fre¬ 
quency  range  of  seismic  signals  (about  0.005  Hz  to  50  Hz)  without  loss  of 
information  due  to  the  band  limiting  properties  of  standard  analog  seismic 
systems  widely  in  use  in  the  1970s.  This  approach  has  led  to  the  estab¬ 
lishment  of  numerous  broad  band  seismological  stations  and  several  broad 
band  networks  and  arrays  worldwide.  The  ORFEUS  data  center  was  founded 
in  1987  and  dedicates  itself  exclusively  to  the  collection  and  distribution  of 
broad  band  digital  records.  Meanwhile,  highly  resolving  analog/digital  con¬ 
verters  have  become  available,  making  the  use  of  broad  band  seismometers 
no  longer  imperative  in  the  pursuit  of  the  goals  mentioned  above.  In  fact, 
broad  band  (BB)  seismograms  can  be  simulated,  making  use  of  band  limited 
digital  registrations  from  two  or  more  different  standard  instruments  at  the 
same  location,  provided  their  combined  pass  bands  cover  the  frequency  range 
^»red  amd  the  resolution  of  analog  to  digital  conversion  is  sufficiently  high. 

The  synthesis  process  involves: 


1.  Restitution  of  true  (albeit  band  limited)  ground  motion  by  deconvolv¬ 
ing  system  response  for  each  instrument. 

2.  Adjustment  of  sample  frequencies. 

3.  Compensation  of  time  shifts. 


4.  Tailoring  individual  pass  bands. 

5.  Arithmetic  summation  of  traces. 

6.  Numerical  differentiation  in  case  a  velocity  or  acceleration  proportional 
BB  seismogram  is  to  be  obtained. 

The  procedure  described  b'  'w  was  developed  and  implemented  for  instru¬ 
ments  installed  at  the  GERESS  regional  seismic  array.  It  is  performed  ex¬ 
clusively  in  the  time  domain  and  uses  forward  recursive  filtering  only,  which 
makes  it  potentially  fast  and  suitable  for  on-line  realtime  applications. 


6.3  An  example 

The  German  Experimental  Seismic  System  (GERESS)  array  project  is  a  co¬ 
operative  research  program,  jointly  undertaken  by  Southern  Methodist  Uni¬ 
versity  in  Dallas,  Texas  and  Ruhr- University  Bochum,  Germany  [Harjes,  1990]. 
The  GERESS  array  is  located  in  the  Bavarian  Forest  area  in  the  southeastern 
part  of  Germany  near  the  border  with  Austria  and  Czechoslovakia.  The  array 
consists  of  25  stations  with  vertical  component  short  period  (1  Hz)  Teledyne 
Geotech  GS-13  type  instruments  sampled  at  40  Hz.  In  addition,  four  of  the 
sites  include  horizontal  component  instruments.  At  the  key  station  of  the 
airay,  GEC2  (A  =  48.84511  N.  <,?  =  13.70156  E)  there  is  a  supplemental 
three-component  set  of  GS-13’s  sampled  at  120  Hz  (high-frequency  element) 
and  a  three-component  set  of  broad-band  seismometers  (Teledyne  Geotech 
BB-13)  sampled  at  10  Hz. 

The  broad  band  au'^  high  frequency  recordings  of  an  event  from  Roma¬ 
nia  were  selected  for  this  demonstration  (see  Figure  6-1.)  This  event  displays 
a  considerable  signal  to  noise  ratio  over  much  of  the  frequency  range  con¬ 
sidered  (see  Figure  6-2.)  It  is  a  **broad  band  event”,  which  could  not  be 
recorded  exhaustively  with  only  one  of  the  two  instruments  involved. 


6.3.1  Restitution  of  true  ground  movement 

Analog  transfer  characteristics  for  the  two  seismometer  systems  GS-IS  and 
BB-13  are  specified  by  the  manufacturer  stating  poles,  zeroes  and  gain  factor 
of  the  transfer  function  in  the  Laplace-domain  (see  table  6-1.) 
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Figure  &-1:  Original  seismograms  recorded  at  GERESS  reference  station 
GSC2,  of  a  broad  band  event  from  Romania  of  July  12,  1991,  NEIS  PDE; 
t  -  10:42:21.3,  v?  =  45.357  N,  A  =  21.052  E,  =  5.0,  A  from  GERESS 
670  km. 
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Figure  6-2:  Smoothed  power  spectrum  of  the  event  from  Romania  (solid 
line)  and  of  noise  immediately  preceding  the  event  (broken  line).  Length  of 
time  window  660  s. 
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Figure  6-3:  Amplitude  (top)  and  phase  (bottom)  characteristics  of  the  Seis¬ 
mometers  GS-13  and  BB-13,  as  specified  by  the  manufacturer  (solid  lines) 
and  approximated  by  poles  and  zeros  below  the  Nyquist  frequency  (broken 
lines) 


Table  6-1;  Gain.  Poles  and  zeroes  of  GERESS  seismometer  transfer  func¬ 
tions.  (For  GS-13:  asstuning  earth  velocity  as  input,  for  BB-13:  assuming 
earth  acceleration  as  input.  Poles  and  zeroes  below  the  Nyquist  frequency 
are  printed  in  bold.  If  these  are  to  be  used  as  an  approximation,  the  gain 
must  be  adjusted  as  indicated  in  brackets. 


Seism.type 

gain  St 

zefoe  poles  f«r 

GS.13 

1.5818E28 

(2.659Sfi-H)l} 

0 

0 

-4.3845E-H)0  +  i  4.1812£-{-00 

-4.3846E-H)0  i  4.1812E-f00 

-1.2218E-M)2  +  i  3.9408E-H)2 

-1.2218E-M)2  -  t  3.9408E-t-02 

-6.0728E-K)2 

-1.3298E-M)3 

-3.5000E-H)7 

-1.7361E-H}8 

60-13 

1.9991E-I-23 

(4.279rE-01) 

-6.1191E~0& 

-3.201SE-M)4 

-1J188E-H}3 

-l.2793E-H)3 

-4.5749wf01  -i-  i  1.3322»4>02 

-4.5749»^mi  -  t  1.2322e-MI2 

-3.0340H-03 

-1.2171e+03 

— 6.0972e+02 

-6.2399e-02 

-l.0388e-f02 

-1.3298«403 

-2.5000H-07 

-1.7381e-M)8 
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Figure  6-4:  Restitution  operators  are  obtained  by  reversing  seismometer 
transfer  functions  (solid  lines.)  They  are  simplified  by  omitting  poles  and 
zeroes  which  have  effects  above  the  respective  Nyquist-frequencies  only,  and 
high  pass  filters  at  T=50s  and  T=300s  are  introduced  to  avoid  instabilities. 
Amplitude  characteristics  of  the  resulting  operators  are  given  by  the  dotted 
broken  lines. 


Assuming  ground  velocity  and  ground  acceleration  as  input  for  GS*13 
and  BB*13^  respectively,  this  leads  to  the  magnification  and  phase  character¬ 
istics  plotted  in  Figure6-3  (solid  lines.)  For  our  purpose,  only  the  intervals 
between  dc  cutoff  filter  (at  3.333  •  10"®  Hz)  and  Nyquist  frequency  (at  5 
Hz  and  60  Hz  respectively)  are  relevant,  thus  the  transfer  functions  can  be 
simplified  by  omitting  all  poles  and  zeros  corresponding  to  frequencies  above 
the  Nyquist  frequency,  retaining  only  those  printed  in  bold  in  table  6-1  and 
adjusting  the  gain  factor  accordingly.  This  manipulation  produces  only  small 
consequences  below  the  Nyquist  frequency  (broken  lines  in  Figure  6-3.) 

Direct  inversion  of  the  transfer  function  making  zeros  to  poles  and  vice 
versa  and  taking  the  inverse  of  the  gain,  results  in  an  operator  capable  of 
exactly  reversing  the  influence  of  the  systems  on  a  signal  (Figure  6-4,  solid 
lines.)  But  as  the  inversion  involves  integration  in  the  case  of  both  systems 
(symbolized  by  poles  at  s  =  0)  these  operators  would  become  unstable  in 
the  presence  of  even  the  smallest  dc  offset.  This  instability  must  be  removed 
by  introducing  two  Butterworth  high-pass  filters  at  7  —  300  s  (BB-13)  and 
7  s:  50  s  (GS-13),  which  can<%l  the  poles  at  zero  (dashed  lines  in  Figure 
6-4.)  The  complete  restitution  oporators  in  the  Liq>lace-domain  are  listed  in 
table  6-2. 

For  digital  application  these  operators  are  transformed  onto  the  z- 
domain  and  applied  as  forward  recursive  filters,  using  sine  tapyers  to  avoid 
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Table  6-2:  Operators  to  restitute  true  ground  displacement  from  GERESS 
high  frequency  channels  (GS-13  sampled  at  120  Hz)  and  broad  band  channels 
(BB-13  sampled  at  10  Hz.)  Operators  are  given  in  gain/ zero/ pole  form  both 
in  the  s-domain  (see  eq.  1 )  and  the  z-domdn  (eq.  13),  as  well  as  in  polynomial 
form  (eq.  14)  in  the  appendices 


Operator 

gain  9s 

zerob  s**k 

poles  5*«*; 

3.7720E-02 

-4.3845E+00  +i  4.1812e+00 
-4.3845E+00  -i  4.18i2e+00 

-1.2566E-01 

-6.2832E-02  +i  1.0883E-01 
-6.2832E-02  -i  1.0883E-01 

gain  g. 

zeros 

poles 

HFINVhpSO 

+9.6354E-0I  +i  3.3602E-02 
+9.6354E-01  -i  3.3602E-02 
-l.OOOOE+00 

+9.9895E-01 

+9.9948E-01  -t-i  9.0642E-04 
+9.9948E-01  -i  9.0642E-04 

1.6284E-04 

numerzitor  coefficients  at 

denominator  coefficients  ii 

+ 1  .OOOOOOOOOOOOOOOE+00 
-9.270707123784236E-0\ 
-9.975412649867432E-01 
-I-9.295294473916804E-01 

+1.000000000000000E-1-00 

-2.997905604993305E+00 

+2.9958134O2657542E-f0O 

-9.979077965171582E-01 

gain  g. 

zeros  s«»t 

poles 

0 

-6.2399E-02 

-6.n91E-05 

-2.0944E-02 

-1.0472E-02  +i  1.8138E-02 
-1.0472E-02  -i  1.8138E-02 

gain  gt 

zeros  zf*t 

poles 

BBINVhpSOO 

+9.937795E-01 

+1.000000E4^)0 

-l.OOOOOOE+00 

-l.OOOOOOE+00 

•f9.999939E-01 

+9.979078E-01 

4'9.989517E-01  +i  1.8U901E~03 
-1-9.989517E-01  -i  1.811901E-03 

S.8474E-12 

numerator  coefficients  at 

denominator  coefficients 

+I.OOOOOOQOOOOOOOOE+00 

-H6.220506692349304E-03 

-1.993779493307651E+00 

-6.220506692349415E-03 

-I-9.937794933076506E-01 

+1  OOOOOOOOOOOOOOOE-t^OO 
-3.99580S091448180E<H)0 
+5.987424068361191E-H)0 
-3.9874328531S6351ENf00 
+9.958138762433976E-01 

transient  effects.  A  procedure  to  effect  the  necessary  transform  is  given  in 
the  Appendix  I,  the  restitution  operators  in  pole/zero  form  in  the  z-domain 
and  as  recursive  filter  coefficients  are  listed  in  table  6-2. 

6.3.2  Adjustment  of  sample  frequencies 

Before  the  two  traces  from  restituted  GS-13  and  BB-13  registrations  can  be 
added,  sample  firequencies  must  be  equalized,  interpolating  1 1  values  between 
every  2  values  in  the  10  Hz  BB-13  trace.  The  nature  of  the  interpolation 
algorithm  employed  is  not  critical,  since  frequencies  falsely  introduced  by 
rough  interpolation  will  be  cut  off  later,  anyhow.  Both  an  elaborate  and  a 
simple  triangular  algorithm  were  tested  and  found  to  perform  equally  well. 


6.3.3  Compensation  of  time  shifts 

At  this  stage  a  possible  time  lag  between  the  two  traces  can  be  found  by 
calculating  the  cross-covariance  function.  In  the  case  of  GEC2,  the  lag  which 
is  already  clearly  visible  on  the  raw  recordings,  was  found  to  be  18  samples 
(corresponding  to  0.13  s.)  The  lag  was  compensated  for  by  shifting  the 
leading  BB-13  trace  with  respect  to  the  trailing  GS-13  trace. 

6.3.4  Tailoring  of  individual  passbands 

If  the  traces  were  added  at  this  stage,  the  middle  frequency  range  -  where 
both  seismometers  perform  well  -  would  be  overrated  as  compared  to  the 
low  frequency  end  below  0,02  Hz  (where  the  GS-13  trace  was  not  properly 
restituted)  and  the  high  frequency  end  above  5  Hz  (where  the  BB-13  trace 
was  merely  interpolated.)  The  solution  is  to  perform  a  weighted  summation 
by  first  applying  a  high-pass  filter  to  the  GS-13  restituted  trace  and  a  low- 
pass  filter  to  the  BB-13  restituted  trace.  These  ‘^tailoring  filters”  are  depicted 
in  Figure  6-5.  Care  is  required  in  their  design  (see  Appendix  II),  because 
they  determine  the  transfer  characteristics  of  the  synthesized  BB  trace  (solid 
line  in  Figure  6-5.)  The  small  imperfection  in  the  amplitude  characteristic 
at  0.2  Hz  is  unavoidable,  when  filtering  is  limited  to  forward  recursion.  The 
tailored  traces  are  shown  in  Figure  6-6. 

6.3.5  Addition  and  numerical  differentiation 

Now  the  weighted  traces  can  be  added  point-wise,  resulting  in  the  desired  BB 
seismogram.  It  is  a  record  of  true  ground  displacement  between  3.333  •  10~^ 
Hz  and  48  Hz  (anti-alias  filter  cut-off.) 

To  obtain  a  velocity-proportional  or  acceleration- proportional  seismo¬ 
gram,  the  BB  trace  may  be  numerically  differentiated  once  or  twice  (Figure 
6-7.) 
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Figure  6-5:  Amplitude  (top)  and  phase  (bottom)  characteristics  imposed  on 
the  synthesized  BB  trace  by  tailoring  filters  (solid  line.)  The  broken  and 
dotted  lines  represent  the  contributions  of  HP  and  LP  filters  used.  Imperfec¬ 
tions  at  the  middle  frequency  are  unavoidable,  they  are  shown  in  linear  scale 
in  fig  0-11. 
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Figure  6-6:  Restituted  seismograms  (interpolated  in  case  of  BB-13),  run 
through  the  "tailoring  filters'',  prior  to  summation. 
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Figure  6-7;  Synthesized  BB  trace.  From  top  to  bottom:  (a)  whole  seismo¬ 
gram  (b)  first  onset  displacement  (c)  velocity  (d)  acceleration 
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Figure  6-8:  Comparison  of  originally  measured  traces  (solid  lines)  and  traces 
recovered  from  BB  synthesis  (broken  lines)  for  verification  purposes. 


6.4  Verification 

For  the  purpose  of  verification  of  the  method,  the  GS*13  original  trace  can 
be  recovered  by  treating  the  BB  trace  as  input  and  simulating  GS-13  system 
response  through  application  of  the  transfer  function  according  to  Appendix 
I.  Figure  6-8  (top)  shows  that  a  reasonable  match  is  achieved  between  the 
measured  and  the  simulated  GS-13  trace.  The  small  differences  visible  in  the 
zoomed  first  few  seconds  can  be  ascribed  to  the  phase  bias  introduced  by  the 
stabilizing  high-pass  filter  at  T  =  50  s  in  the  restitution  phase. 

The  same  should  theoretically  be  possible  for  the  BB-13  original  trace, 
but  here  the  anti-alias  filter  (3  dB  point  at  4.5  Hz)  of  the  digital  recording 
system  must  be  modelled  carefully  and  the  number  of  data  points  decimated 
by  a  factor  of  12.  This  introduces  a  certain  amount  of  arbitraryness  (at  which 
point  to  start  decimation)  and  the  coincidence  is  not  as  good  (see  Figure  6-8, 
bottom.)  The  match  could  be  improved,  however,  by  low-pass  filtering  both 
the  original  and  the  recovered  trace. 

One  can  therefore  be  reasonably  sure  that  after  application  of  the  above 
procedure  one  has  a  faithful  BB  record  of  true  ground  movement,  which  can 
be  used  in  further  studies  requiring  broad  band  data  (e.g.  moment  tensor 
inversion),  or  can  be  filtered  to  match  the  system  response  of  seismograph 
systems  serving  as  international  reference. 
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Figure  6-9:  Tailoring  filters  which  were  applied  in  the  frequency  domain 
by  Harvey  &  Choy  (1982].  While  adding  up  to  unity  in  the  frequency  do¬ 
main  (dashed  lines),  their  discontinuous  derivatives  should  produce  sinc^(t)- 
ripples. 
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6.5  Discussion 

Problems  may  arise,  if  the  signal  to  noise  ratio  drops  below  unity  in  the  low 
frequency  range.  In  this  case  the  long  period  noise  would  be  amplified  in 
the  BB  trace,  completely  overwhelming  the  short  period  signal.  Rather  than 
redesigning  the  restitution  filters  with  an  increased  high-pass  cut-off  point  it 
is  recommended  to  prefilter  the  original  traces  with  a  high-pass  filter  before 
the  application  of  the  BB  synthesis.  Alternatively  the  BB  trace  could  be 
high-pass  filtered. 

The  small  phase  bias  introduced  by  the  restitution  filters  could  be 
avoided,  if  two-sided  recursive  algorithms  (described  by  Marschall  [1984]) 
were  used  for  the  implementation  of  all  filters.  This  would,  however,  pre- 
-  dude  the  on-line  real-time  application  of  the  synthesis. 

The  same  applies  for  synthesis  in  the  frequency  domain,  following  Har¬ 
vey  &  Choy  [1982].  Moreover,  the  overlapping  frequency  filters  employed  by 
these  authors  (see  Figure  6-9.)  have  discontinuous  deri  vatives,  which  should 
produce  sinc^(t)-ripple8  in  the  time  domain. 

The  described  method  can  also  be  applied  to  narrow  band  seismograph 
systems  like  WWSSN-lp  and  WWSSN-sp,  as  long  as  the  resolution  of  the 
aoalog/digital  converter  is  suffidently  high  to  raise  the  magnification  to  unity 
in  the  overliq>ping  region  while  preserving  an  acceptable  signal/noise  ratio. 
As  can  be  seen  on  Figure  6-10  an  A/D-converter  with  24  bit  resolution,  as 
used  at  GERESS,  would  suffice  in  this  case. 
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Figure  6-10:  Amplitude  curves  of  WWSSN-Ip  (solid  line)  and  WWSSN- 
sp  (broken  line)  system  response  to  ground  velocity.  [Static  gain  factors 
were  adjusted  to  1,000  and  200.000  for  Ip  and  sp  respectively]  In  order  to 
synthesize  a  BB  (10~^ .  ..20Hz]  trace,  the  signal  level  must  be  raised  by  up 
to  2.5  decades.  In  case  of  a  24-bit  A/D-converter  as  used  at  GERESS,  this 
would  leave  at  least  95  dB  resolution  over  the  whole  frequency  range. 
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APPENDIX  I 


Let  H{s)  be  a  stable  transfer  operator  in  the  ^-domain,  s  =  iu>  being  the 
variable  of  the  Laplace  transform.  It  is  usually  given  in  (or  can  be  cast  into) 
pole/zero/gain  forr^ 


ffla\  —  n  „ 


n  (a  -  sWfc) 

ksO 

N 

n  (s  -  sM() 

IsO 


with  gain  g^,  zeroes  sWfc  and  poles 

In  order  to  be  applicable  to  time  sampled  data,  it  must  first  be  approx¬ 
imated  in  the  r-domain,  z  =  being  the  variable  of  the  Z-transform  with 
sampling  interval  T.  Several  transformations  were  tried  (pulse-invariant, 
step-invariant  and  ramp-invariant,  see  (Schufiler,  1964]  for  a  review)  but  the 
well-known  bilinear  z-transformation  [Golden,  1964] 

2  2-1 

»=r-7?T 


together  with  frequency  prewarping 


-  2  /T 


(I-) 


was  found  to  perform  best. 

Single  poles  and  zeros  are  simply  prewarped  according  to 

2  /T  „  \ 


5«voo»  =  _  .  tan 
T 


(l-H- 


Complex  conjugate  poles  and  zeros  can  be  imagined  to  arise  as  solutions  of 
quadratic  equations  of  the  form 

s’ +  2hun)S  +  tuo*  =  0  (5) 

which  describe  damped  harmonic  oscillation  in  the  s-domain  with  free  circu¬ 
lar  frequency  ujn  and  damping  constant  h  <  1.  Their  solutions  are 


Introducing 


Si/2  =  •  [a  ± » ( 1  -  h’)*'’] 

X  :=  3S(Si)a:  — WqA 
y  ;ss  9(5,)  ss  — u^oV^l  —  h? 


and  solving  for  uq  and  h  we  get 

UAo  =s  y/x^  -b  y*  *|«j| 

X  |»(st)| 
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h  is  left  unchanged,  while  wo  is  prewarped  to  wq  according  to  (3),  such  that 
the  prewarped  conjugate  pair  is 

[a  T  t(l  “  .  (9) 


In  all  the  bracketed  expressions  in  eq.  1,  the  s<‘^>  is  replaced  by  sw*' 
fronm  eq.  4  or  eq.  9,  while  s  is  replaced  by  y  •  ~  from  eq.  2.  Thus  each 

and  by  rearranging 


2,  _  jw-).  /  2  +  s<f>^h-T\ 

z  —  I  V  2  —  swoo, .  T ) 


Identifying 


2  +  J«Vbe);  .  T 
2  —  -  7* 


T(<V80)j. 


and 


(11) 


we  get 


3 

the  desired  operator  in  the  ;-domain. 

For  each  pole  and  zero  in  the  s-domain  a  corresponding  pole  or  zero  in 
the  2-domain  can  be  computed  according  to  eq.  11  and  the  gain  is  adjusted 
according  to  eq.  12.  For  every  excessive  pole  {N  >  M)  in  the  s-domain, 
one  zero  at  r  =  —  1  is  added  in  the  z-domain.  In  case  of  excessive  zeros, 
(z  +  1)^"*^  would  produce  one  or  more  poles  on  the  unit  circle,  thereby 
destabilizing  the  filter.  This  shows,  that  the  transform  is  not  suitable  in 
case  of  excessive  zeroes.  Another  condition  derived  from  eq.  3  is  that  no 
characteristic  frequency  associated  with  a  pole  or  zero  via  eq.  7  may  be  higher 
than  the  Nyquist  frequency  ^  because  of  the  periodicity  of  the  tangent 
function. 

Even  though  these  conditions  are  met,  the  accuracy  to  which  H{z) 
i4>proximates  if(s)  varies,  being  usually  poor  close  to  the  Nyquist  frequency. 
It  can  be  checked  by  evaluating  H{z)  along  the  unit  circle  and  comparing  it 
to  ir(s)  evaluated  along  the  imaginary  axis. 
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Once  eq.  13  has  been  obtained,  it  can  be  brought  into  polynomial  form 


=  -  (14) 

E 

1=0 

by  multiplication  and  is  implemented  as  a  forward  recursive  filter  following 
[Seidl,  1980,  eqs.  10  through  16)  or  (Ferber,  1989,  eqs.  26  and  48  through  51] 

APPENDIX  II 

Overlapping  low  pass  and  high  pass  Butterworth  filters  were  employed  as 
tailoring  filters  because  of  their  smoothness  both  in  the  time  and  frequency 
domain.  Their  amplitude  characteristics  Mip{u))  and  Mppiui)  are.  given 
[Kulhanek,  1967,  pp  87ff]  by 


Mlp{u) 


uijj*  +  (  — l)"ti; 


in 


M^p(u}) 


(- 1 


(15) 


They  are  solely  determined  by  orders  n  and  comer  frequencies  and  afp, 
which  shall  be  spaced  equally  (on  a  logarithmic  frequency  scale)  around  a 
middle  frequency  ujm- 


=  u/M  •  a  .  (16) 


The  condition  that. 


/ 

Mi,(u  =  wm)  +  =  o)m)  =  1 

and,  because  of  symmetry, 

=»  —  Mh{^m)  =  0.5  :=  a 

leads  to  the  expression  for  a 


a  = 


(17) 


With  set  to  2jr  •0-2H2  and  the  ord^  n  =  4  chosen  for  coincidence  of 
phase  characteristics  of  both  LP  and  HP,  comer  frequencies  can  be  calculated 
according  to  eq.  16.  Unfortunately,  the  resulting  amplitude  characteristics 
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amp.  (linear) 


Figure  0-11:  Close-up  of  the  crossover  point  of  tailoring  filters  in  linear 
scale.  Ideally  the  sum  of  the  filter  amplitudes  should  be  unity  everywhere. 
The  condition  in  eq.  17  yields  unity  at  0.2  Hz  (left),  while  the  condition  in 
eq.  18  results  in  an  amplification  which  preserves  overall  gain  (right.) 


frequency  [Hz], linear 


(Figure  0-11,  left)  do  not  add  up  to  unity  along  the  whole  frequency  axis.  If 
over2dl  gain  is  to  be  preserved,  the  expression 

7  f 

J  [A/ti>(w,a)  -h  (A/ff/>(uJ,a)  -  1]  dw  *  0  ,  (18) 

0 

evaluated  by  numerical  integration  ana  a  zero-finding  algorithm,  leads  to  the 
determination  of  an  optimum  value  of  a  (n  4,  wl  ==  2x  •  0.2  Hz)  =  1.19504. 
The  corresponding  amplitude  characteristic  is  also  plotted  in  Figure  O-Il 
(right.)  In  the  example  considered,  the  resulting  differences  in  the  BB  trace 
were  below  3%. 

Poles  and  zeros  of  the  filters  were  calculated  following  [Kulhanek,  1967, 
eqs.  3.16  and  3.17].  Transfer  to  the  z-domain  and  implementation  as  forward 
recursive  filters  were  done  as  described  in  Appendix  I. 


/  thank  Prof.  Harjes,  who  suggested  the  topic,  and  my  colleagues  for 
providing  a  creative  working  environment. 
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7.  QUANTIFICATION  OF  TIME  -  INDE¬ 
PENDENT  SPECTRAL  STRUCTURES 
IN  SONOGRAMS 

Jan  Wiister 


7.1  Introduction 

Near  regional  seismograms  of  conventional  industrial  explosions  (i.e.  mine 
and  quarry  blasts)  show  a  pecularity  termed  “spectral  scalloping''  or  “time- 
independent  spectral  structures",  which  is  absent  from  earthquake  seismo¬ 
grams.  It  consists  of  local  spectral  peaks  and  troughs  that  can  be  traced  along 
the  whole  duration  of  the  seismogram.  Sonogram  analysis  and  subsequent 
smoothing  and  l-bit-compression  bring  out  these  structures  as  horizontal 
batnds  (Hedlin  et  al.,  1989).  The  phenomenon  is  observed,  whenever  indus¬ 
trial  explosions  are  studied  using  high  frequencies,  it  is  generally  explained 
with  the  multiple  nature  of  commercial  blasts,  the  so-called  ripple-firing  tech¬ 
nique  (peaks)  (e.g.  Smith,  1989)  and  by  the  finite  size  and  duration  of  the 
source  (troughs)  (e.g.  Chapman  et  al.,  1992.) 

The  presence  of  spectral  structures  has  been  used  as  a  discriminant 
between  earthquakes  and  chemical  explosions  in  discrimination  studies  in 
Scandinavia  using  NORESS  data  (Baumgard  k.  Ziegler,  1988,  Dysart  k  Pulli. 
1990)  and  also  in  the  Vogtland  area  in  Central  Europe  using  GERESS  data 
(Wuster,  1993.)  In  Figures  7-1  and  7-2  two  typical  examples  from  the  earth¬ 
quake  and  explosion  populations  of  the  latter  study  are  compared.  Both 
earthquakes  and  explosions  display  spectral  complexity  in  the  sonograms. 
However,  only  in  the  case  of  the  blasts  can  the  local  peaks  be  traced  along 
the  whole  signal  portion  of  the  sonogram  between  the  P-onset  at  t=10  sec 
and  the  coda  at  approximately  tssSO  sec,  as  revealed  by  the  1- bit-compressed 
digital  matrices. 

7.2  Why  quantify? 

The  presence  or  absence  of  time-spanning  spectral  structures  is  easily  noted 
in  Figures  7-1  and  7-2,  but  in  some  cases,  particularly  with  poor  signal/noise 
ratio,  the  distinction  is  not  as  clear-cut.  Also,  in  a  concept  in  which  several 
discriminants  are  combined  to  produce  a  final  verdict,  quantitative  values 
must  be  attached  to  each  discriminant.  Moreover,  if  the  origin  of  the  spectral 
structures  is  attributed  to  ripple-firing,  characteristic  ‘wavelengths'  in  the 
frequency  domain  are  indicative  of  the  delay  between  successive  explosions 
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Figure  7-1:  Seismogram,  sonogram  and  1-bit-compressed  digital  matrix  of  a 
typical  Vogtland  earthquake 
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Figure  7-2:  Seismogram,  sonogram  and  1-bit-compressed  digital  matrix  of  a 
typical  Vogtland  explosion 
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and  may  aid  in  the  identification  of  specific  quarries  or  mines.  Therefore  a 
quantitative  evaluation  of  the  spectral  structures  is  desirable.  To  this  end 
the  following  techniques  were  applied: 

•  cepstral  analysis 

•  harmonic  product  spectrum 

•  horizontal  sums. 


7.3  Cepstral  analysis 


Let  f{t)  be  a  time  series  and  |F{u;)j  a  crude  estimation  of  its  amplitude- 
spectrum  derived  from  the  Fourier  transform  =  F {(*>)•>  then  the  real 

cepstrum  /(f)  is  defined  according  to  Bath  (1974)  as  the  amplitude-spectrum 
of  the  real  logarithm  of  |F(u;)|: 


fit)  = 


+O0 

I 


log(iF(u;)l)  e-'^du, 


(1) 


Its  merit  as  opposed  to  a  simple  spectrum  of  the  spectrum  lies  in  the  fact, 
that  the  logarithm  seperates  the  different  factors,  if  the  spectrum  can  be 
written  as  a  product 

Fiu)  =  Giu>)  •  H{u)  .  (2) 

This  is  regularly  the  case  with  multiple  sources,  as  shown  in  the  synthetic 
example  of  Figure  7-3,  Five  successive  2nd-order  Kupper-pulses  of  0.1  sec 
duration  each  are  applied  with  delays  of  0.2  sec  (Fig.  7-3a.)  The  result¬ 
ing  spectrum  |/’(4*?)|  (Fig.  7-3b)  is  the  spectrum  of  the  single  Kupper-pulse 
with  its  spectral  maximiun  at  10  Hz  (|G(u;)|,  Fig.  7-3c)  multiplied  by  the 
spectrum  of  five  successive  spikes  (|/f(u>)|,  Fig.  7-3d.)  The  spectrum  of  the 
spectrum  (Fig.  7-3e)  is  -  according  to  the  convolution  theorem  -  the  con¬ 
volution  g{t)  *  h{t).  The  peak  at  0.2  sec  in  Figure  7-3e  is  accompanied  by 
two  side  maxima  which  would  tend  to  broaden  and  obscure  it  in  the  presence 
of  noise,  neither  could  the  second  main  frequency  component  at  0.1  sec  be 
determined  in  an  easy  way.  In  the  cepstrum  (Fig.  7-3f),  on  the  other  hand, 
both  spectral  factors  are  clearly  seperated  into  distinct  spectral  peaks. 

This  property  of  the  cepstrum  is  used  for  vocal  pitch  detection  in  acous¬ 
tics  (Noll,  1963)  and  a  corresponding  apparatus  for  use  in  telephony  has  even 
been  patented  (Schroeder,  1965.)  The  discrimination  between  single  events 
and  ripple-fired  explosions  in  seismology  poses  an  analogous  problem  to  the 
discrimination  between  consonants  and  vowels  in  speech  processing.  In  a 
voiced  portion  of  speech,  the  vocal  tract  is  excited  by  periodic  puffs  of  air 
produced  by  the  glottis,  just  as  in  a  ripple-fired  explosion  the  wave-path  to 
the  receiver  is  excited  by  periodic  explosions.  The  analogy  suggests  that  the 
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Figure  7-3:  e,f 


cepstrum  should  be  well  suited  to  extract  the  dominaat  quehrency  (swave- 
length  in  the  frequency  domain),  thus  quantifying  the  effect  of  the  ripples. 

Unfortunately,  the  logarithm  is  very  susceptive  to  noise,  which  is  much 
more  a  problem  in  sdsmology  than  in  speech  processing.  This  is  shown  in 
Figure  7-4  with  the  same  synthetic  example  adding  only  10%  of  normally 
dirtributed  random  noise.  While  the  periodicity  of  the  spectrum  is  still  very 
clear,  the  corresponding  cepstral  peak  is  already  quite  feeble,  and  with  33% 
noise  added  it  is  no  longer  recognizable. 

In  cases  of  high  signal/noise  ratio,  the  cepstrum  can  indeed  help  to 
quantify  frequency  structures  in  sonograms,  as  is  shown  with  a  synthetic  ex¬ 
ample  in  Figure  7-5.  The  time  series  RIPP120  is  made  to  resemble  the  data 
traces  used  in  the  discrimination  study.  It  is  made  up  of  Kupper  pulses  Of 
0.083  sec  duration  each,  at  0.167  sec  interval  from  f  s  10 ...  45  sec;  1%  nor¬ 
mally  distributed  random  noise  is  superposed.  Sonograms  are  calculated  us¬ 
ing  exactly  the  same  procedure  as  with  actual  data.  In  the  1-bit-compressed 
matrix  (Fig.  7-5a)  a  clear  banding  structure  results.  The  spectral  period  read 
from  the  matrix  is  6  Ebs.  This  corresponds  to  a  quefrency  of  l/6Hz  »  0.16$ . . . 
sec  brottj^t  out  by  the  dominant  cepstral  peak  in  the  wiggle  trace  plot  of 
Figure  7-5b. 

Turning  now  to  real  data,  the  cepstral  wiggle  trace  plot  for  the  event 
VS014  (Fig.  7-6a)  does  not  allow  a  dominant  quefrency  to  be  pidked,  al¬ 
though  the  onxesponding  banding  structure  in  Figure  7-lb  is  also  very  clear. 


Ill 


Figure  7-4:  a,b,c,d:  Synthetic  example  with  noise 
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Figure  7-5:  Synthetic  data  with  1%  noise:  a)  1-bit-compressed  digital  matrix, 
b)  cepstrum  matrix 


Against  these  odds,  Baumgard  k.  Ziegler  (1988)  have  successfully  em¬ 
ployed  cepstral  peaks  for  discrimination,  using  spectral  stacks  of  all  NORESS 
channels  for  smoothing  and  noise  reduction.  This  is  not  an  option  here,  if 
frequencies  above  16  Hz  are  to  be  considered,  because  the  normal  GERESS 
channels  are  sampled  at  40  Hz  and  anti-alias  filtered  at  16  Hz.  Just  the 
hf-element  at  array  station  GEC2  is  sampled  at  120  Hz  (anti-alias  filter  at 
48  Hz.) 


7.4  Harmonic  product  spectrum 

Another  method  for  the  measurement  of  the  pitch  or  fundamental  frequency 
of  a  signal  is  the  harmonic  product  spectrum  (HPS).  It  is  defined  by  Schroeder 
(1986)  as  the  function  obtained  by  adding  the  logarithm  of  all  spectral  am¬ 
plitudes,  which  are  integer  multiples  of  the  fundamental  frequency,  i.e.  the 
harmonic  overtones: 

,V  N 

HPS(u;)  Y.  =  log  n  •  (3) 

nsl  nsl 

The  HPS  should  give  a  clear  peak  at  the  fundamental  frequency,  even  if  this 
frequency  and/or  some  of  the  overtones  are  missing,  provided  the  overtones 
present  are  harmonic.  The  S3mthetic  example  in  Figure.  7-3  has  already 
shown  that  periodic  excitation  results  in  harmonic  overtones,  and  the  funda¬ 
mental  frequency  of  5  Hz  is  indeed  extracted  well  by  the  HPS  in  Figure  7-7. 
even  in  the  presence  of  noise. 

Application  of  the  HPS  to  the  sonogram  matrix  for  the  synthetic  ex¬ 
ample  RIPP120  yields  a  clear  peak  at  the  fundamental  frequency  of  6  Hz 
(Fig.  7-7c).  However,  application  of  the  method  to  real  data  does  not  yield 
a  comparable  peak  (Fig.  7-6b) 

7.5  Horizontal  sums 

At  this  point  it  was  decided  to  start  from  the  phenomenology:  If  spectral 
bands  are  observed,  surely  they  can  be  quantitatively  extracted.  Other  than 
in  the  applications  of  pattern  recognition  methods  to  the  detection  of  seis¬ 
mic  horizons  (e.g.  Lu  &  Cheng,  1990),  we  expect  bands  in  the  sonograms, 
corresponding  to  time^independent  spectral  peaks,  to  be  strictly  parallel  and 
horizcmtal. 

In  <Mrder  to  extract  such  structures,  horizontal  sums  are  calculated  ac- 
cross  the  l-bit-compressed  digital  sonogram  matrices,  along  the  time  axis. 
White  spots  (local  spectral  highs)  are  assigned  a  value  of  1,  black  spots 
(local  spectral  lows)  are  assigned  a  value  of  0.  In  order  to  avoid  the  misin- 
terpretatitm  of  stationary  spectral  peaks  of  the  noise  as  spectral  bands  of  the 
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Figure  7-6:  Cepstrum(a)  and  Harmonic  Product  Spectrum(b)  matrices  for 
a  clear  Vogtland  explosion 
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Figure  7-7:  Harmonic  Product  Spectrum  (HPS)  of  synthetic  data 
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Figure  7-8:  Horizontal  sums  accross  the  l-bit-compressed  digital  matrices. 
(Solid  lines:  over  signal  sections,  broken  lines:  over  noise  sections.)  Upper 
trace:  Vogtland  earthquake,  lower  trace:  Vogtiand  explosion 


signal,  seperate  sums  are  calculated  for  the  signal  section  (t  ~  10. . . GOsec) 
and  the  noise  sections  (t  =  0 ...  10  and  60 . . .  100  sec)  of  each  matrix.  Two 
examples  are  shown  in  Figure  7-8.  A  band  is  declared,  where  the  signal 
sum  surpasses  the  75%  threshold  and  the  noise-sum  does  not.  A  tool  was 
developed  which  enables  the  observer  to  interactively  determine  the  middle 
frequency  of  these  bands  (Fig.  7-9.)  In  the  Vogtiand  discrimination  study 
mentioned  above,  the  number  of  bands  thus  determined  discriminated  effec¬ 
tively  between  a  population  of  commercial  blasts  and  earthquakes,  as  shown 
by  the  histogram  in  Figure  7-10. 


7.6  Determination  of  blasting  delays 

A  study  on  Czech  quarries  has  given  reason  to  expect,  that  most  of  the 
explosions  used  in  the  Vogtiand  discrimination  study  ori^nate  in  the  two 
coal  strip  mines  "Jiri  Vintirov”  and  *‘Nove  Sedlo**,  within  2  km  of  each  other, 
close  to  the  town  of  Karlovy  Vary.  Speculating  that  blasting  techniques  may 
therdcffe  be  similar  in  many  cases,  the  detected  spectral  bands  were  stacked 
in  spectral  bins  of  1  Hz  width.  The  resulting  histogram  is  shown  in  Figure  7- 
11,  the  positicMDS  of  statkmary  noise-bands  are  alto  indicated. 

.  Peaks  of  4  or  more  are  obtained  from  the  iustof^am  at  frequency  bins 
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Figure  7-9:  A  tool  to  interactively  pick  the  frequency  bands 


Figure  7-10;  The  number  of  frequency  bands  observed  is  a  powerful  discrim¬ 
inant  between  earthquakes  and  blasts  in  the  Vogtland  region 
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Figure  7-11:  Stacked  spectral  bauds  for  Vogtland  explosions 
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of  2,  10,  18,  25  and  31  Hz.  The  distance  between  successive  peaks  is  fairly 
regular  (6«8  Hz,  on  the  average  7.5  Hz)  There  is,  however,  no  fundamental 
frequency,  of  which  the  detected  bands  are  integer  multiples.  This  is  to  say 
we  are  looking  at  unharmonic  overtones.  The  spectral  periods  observed  cor¬ 
respond  to  quefrendes  between  125  and  166  msec  (138  msec  on  the  average.) 
While  this  is  clearly  above  typical  inter-hole  ddays,  it  is  of  the  order  of  typ¬ 
ical  inter-row  delays  reported  by  Chapman  et  al.  (1991)  from  Appalachian 
explosions.  Until  now,  no  contact  has  been  established  to  the  Czech  blasting 
engineers  in  order  to  compare  actual  blasting  techniques  to  these  results. 

7.7  Discussion 

In  theory,  ripple  firing  shotild  give  rise  to  spectral  reinforcements  at  harmonic 
overtones  to  the  fundamental  frequency,  which  corresponds  to  the  inter-shot 
delay.  A  number  of  authors  has  noted  that  this  phenom<mon  is  very  suscep¬ 
tible  to  disturbances  almost  inevitable  in  blasting  practice.  Chiq>man  et  al. 
(1991)  have  shown  that  the  interference  of  delays  between  mdividual  holes 
and  the  delays  between  rows  of  h<^es  already  introduces  substantial  noise 
into  the  spectrum.  Smith  (1989)  argues  that  the  main  effect  of  time  scatter 
in  blasting  c^  would  be  to  suppress  higher-<xder  harmonics,  leaving  lower- 
order  harmonics  intact.  But  Reamer  et  al.  (1992)  obtain  severe  suppression 
and  drift  of  maxima  already  at  a  rime  scatter  of  only  2%.  At  6%  time  scat- 


ter,  a  realistic  figure  according  to  tbese  authors,  even  the  first  reinforcement 
is  shifted  from  its  theoretical  value.  In  addition  the  effects  of  scatter  in  the 
charge  sizes  must  be  considered. 

The  authors  (Reamer  et  al.  1992)  were  able  to  observe  the  harmonic 
spectra  predicted  by  theory  within  100  m  of  small  ripple-fired  explosions,  us¬ 
ing  a  precise  electronic  firing  system.  At  regional  distances  from  complex  ex¬ 
plosions  fired  with  conventional  (chemical)  delays,  only  unharmonic  and  noisy 
spectral  reinforcements  can  be  observed  in  the  form  of  time-independent  spec¬ 
tral  structures.  Under  these  circumstances,  cepstral  and  harmonic  methods 
may  be  inadequate  to  quantify  the  phenomenon.  The  calculation  of  horizon¬ 
tal  sums  accross  1 -bit-compressed  sonogram  matrices  explicitly  makes  use  of 
the  time- independent  nature  of  the  structures  and  is  well  suited  to  quantify 
the  observed  spectral  banding. 


The  idea  to  try  horizontal  sums  emerged  in  a  discussion  xeith  Johannes 
Schweitzer,  which  is  gratefully  acknowledged. 
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